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Abstract. We study the linearization of the Dirichlet-to-Neumann map for Poincare-Einstein 
metrics in even dimensions on an arbitrary compact manifold with boundary. By fixing a suitable 
gauge, we make the linearized Einstein equation elliptic. In this gauge the linearization of the 
Dirichlet-to-Neumann map appears as the scattering matrix for an elliptic operator of 0-type, 
modified by some differential operators. We study the scattering matrix by using the 0-calculus 
and generalize a result of Graham for the case of the standard hyperbolic metric on a ball. 



1. Introduction 

Let M be the interior of a compact {n + l)-dimensional smooth manifold M with non-empty 
boundary dM. Let x be a fixed smooth defining function for 9Af, i.e. x S C°°(M), x > on M, 
x~^(0) = dM and dx\gj;j ^ 0. Let m e Nq and < a < 1. We say a metric g+ on M is C™'" (resp. 
C™, C°°) conformally compact if x^g+ extends continuously to a (7™'" (resp. C™, C°°) metric on 
M , denoted by For a C° metric 7 on dM , the conformal class of 7 is the set of all metircs on 
dM which is conformally equivalent to 7, denoted by [7], i.e. 

[7] = {e2"7 : u e C°{dM)}. 

If (7+ is c™'" conformally compact, [x'^g+\rpgjj] is said to be the conformal infinity of 

Suppose that n > 3 and odd. For < a < 1, a C^'" Poincare-Einstein metric is a C^'" 
conformally compact metric on M with conformal infinity [7] which satisfies the condition 

Ric{g+) = — ?^.9+ on Af, 
x^g+WdM e [7] on i9M. 

If .g+ is a Poincare-Einstein metric with conformal infinity [7] , then for a choice of representative 
9o G there exists a defining function p such that for some e > and an identification of a 
neighborhood of dM with [0, e) x dM , g+ takes the geodesic normal form 

(2) g+^ p-^idp" +g{p)) 

where g{p) is a 1-parameter family of metrics on dM satisfying g{Q) — go. Here p is called the 
geodesic defining function for (g^,go). If go is C°° , a boundary regularity theorem shows that 
g{p) e C°°([0,e) X dM;S'^T*dM) (See [Al][A4][CDLS][H]). Moreover the Taylor expansion of g{p) 
is even to order n, i.e. 

(3) g{p) = go + p'g2 + ... + p"-\g„-i + p".g„ + 0(p"+^) 

where gi G C°° {dM , S'^T* dM) for i = 2,4, l,n, and (?„ is transverse-traceless on {dM,go), 

i.e. 

(4) trgog„ = 0, Sg„gn = 0, 
where in any local coordinates {U;y^, y") of dM, 

trgn5„ = [go]'''[gn]kl, Sgf,9n = [goY'^igrillj-kdy'- . 

Notice that for i = 2, 4, n — 1, gi is locally determined by the boundary metric go; however, (?„ is 
locally undetermined subject only to the two conditions in The terms go and g„ play the role 
of Dirichlet data and Neumann data respectively in this problem. 

The Dirichlet-to-Neumann relation JV is defined as the set of pairs (0, 1/') such that there exists a 
Poincare-Einstein metric gj^ on M with geodesic normal form and Taylor expansion ([3]) such that 
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go — 4>, gn — ip- ^ satisfies the following equivariance properties with respect to diffeomorphisms 
and conformal changes: 

(5) (^,v^)e^^(<i>V,**^)e^, $=¥|aM, v¥eDiff(M), 

(6) (0, V) e ^ (e^"0, e(^-")"V') Vm G (dM). 

If for a fixed boundary metric 0, the Poincare-Einstein metric 5+, with go = (j) in the geodesic 
normal form, is uniquely determined up to diffeomorphism restricting to identity on the boundary, 
then the relation ^ defines the Dirichlet-to- Neumann map on the boundary denoted hy M : (j) ^ ^. 
However, the global existence and uniqueness of Poincare-Einstein metrics fails for a general compact 
manifold M. 

In this paper, we study the linearization of the Dirichlet-to- Neumann map under the assumption 
that it is well defined. For this goal, we do not need the full strength of uniqueness, the global 
uniqueness, of Poincare-Einstein metric, but only need to assume the uniqueness of the solution to 
the linearized equation at a fixed Poincare-Einstein metric (7+, 

(7) RiCg^h + nh = 

up to the deformation generated by a vector field vanishing on the boundary. See more details in 
Section [51 

For < a < 1, let PE'^-°'{M) C x-'^C'^-'^iM] S^T*lA) be the subspace of C™'" Poincare- 
Einsten metrics on M, which can be given the (7™'" topology on M via a fixed compactification 
as in ([T]). Obviously, this topology is independent of the choice of smooth defining function. Let 
Mei°°(9M) denote the space of C°° metrics on dM. Clearly, the tangent space of Met°°{dM) is 
C°°{dM; S^T*dM). Hence, our Dirichlet-to-Neumann map and its linearization, if well defined, are 

N : Met°°{dM) — > C°°{dM; S^T*dM), 

dJ\f : C°°{dM; S'^T*dM) — > C°°(dM] S^T*dM). 

To set up the problem, let go G Met°°{dM) be a representative of the conformal infinity of some 
C°° conformal compact Poincare-Einstein metric g+. Assume that there exists a neighborhood W 
of go, W C Met°°{dAI), such that for each go S W, there exists a unique Poincae-Einstein metric 
g+ e PE^'", up to diffeomorphisms restricting to the identity on dM, such that x'^g+\j,gjj g [^0]. 
Then the Dirichlet-to-Neumann map is well defined on W. Let 'tig„ho' be the trace-free part of a 
symmetric 2-tensor ho with respect to go- Our main theorem is the following: 

Theorem 1.1. The linearization of the Dirichlet-to-Neumann map at go, dM , is a pseudo-differential 
operator of order n with principal symbol 

a„(dAA) - 2-"^^a„(A|„(tf3„ -tf,„5;(5,„tfg„,5;)-i<5,„tfgj) 

where ho G {dM; S'^T*dM). Moreover, if 

S = {/lo e C°^{dM; .y'^T*dM) : tig^ho ^ 0, Sg„ho = 0} 

is the subspace of trace-free and divergence-free symmetric 2-tensor fields, then dA/'ls is elliptic with 
principal symbol 

a„(d^A|s) = 2-"^^a„(A,=J. 

Our work is inspired by the short paper [Grl] of Graham, who studied the Dirichlet-to-Neumann 
map in the case M = B"+^ = {y = (2/°, y") e : ELo l^i^ < !}■ Choosing a defining 

function p = i(l — |yp), the hyperbolic metric <?+ = X)r=o('^y')^ ^ Poincare-Einstein metric 
with prescribed infinity P^5+|ts" = 9o, where go — is the standard sphere metric. Based on his 
earlier joint work [GL] with Lee, the Poincare-Einstein metric uniquely exists with prescribed infinity 
sufficiently closed to go, up to the diffeomorphisms restricting to identity on the sphere. Therefore, 
the Dirichlet-to-Neumann map is well defined in a neighborhood of go in C'^'"(S'"; S'^T*S'") for 
k > n. Graham gave a precise formula for the linearization of the Dirichlet-to-Neumann map at 
go = \d9\'^, which is a pseudo-differential operator and especially, is elliptic on the trace- free and 
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divergence-free sections of (0,2)-tensors with respect to Qq. The generaUzation here of Graham's 
result does not give an exact formula, but only the principal symbol for this operator. See Section 

m 

The outline of this paper is as foUwos. In Section 2, we set up the problem and reduce it to the 
study of two 0-type differential operators. In Section 3, we introduce the 0-calculus due to Mazzeo 
and Melrose, which is used to deal with the geometric operators associated to some asymptotically 
hyperbolic metric in the interior of a compact manifold and is the main technique used here since 
a Poincare-Einstein metric is clearly asymptotically hyperbolic. In Section 4, we choose a suitable 
gauge, the Bianchi gauge, to make the linearized Einstein equation elliptic in the 0-calculus and 
write out the linearized relation of a family of Poincare-Einstein metrics in this gauge and in the 
geodesice normal form. In Section 5, using the method introduced in [MM], [Ma], [GZ] and [JS], we 
construct the parametrix, Poisson kernel and scattering matrix for the linearied Einstein equation in 
the chosen gauge. In Section 6, we prove the main theorem and conclude that the scattering matrix, 
modified by some differential operators, provides the linearization of the Dirichlet-to-Neumann map. 
Finally, in Section 7, we also deal with the N point problem for > 2, to give a better description 
of the Dirichlet-to-Neumann map following the constructions in Sections 2-6. 

I would like to thank Richard Melrose for his help on this paper. The many discussions we had on 
these topics were indispensable and without his patience and guidance, this paper would not exist. 

2. Geometric Setting and Linearization 

Let X be the fixed smooth defining function as in Section [T] Suppose the dimension of DM, n, 
is > 3 and odd. For any q G dM, let {U; y^, y") be the local coordinate patch around q in dM. 
Then there exists e > 0, such that ([0,e) x U;x,y^, ■■•,?/") is a local coordinate patch around q in 
M. For simplicity, denote hy y° = x and {y^joKi^Kn = {x,y) with y = (?/\ 

Suppose that for some < a < 1, (7+ is a C^'" conformally compact Poincare-Einstein metric 
on Af, for which there is a smooth representative, go, in the conformal infinity [x^g+l^g-j^]. Let 
p = p{x, y) be the geodesic defining function for (17+, go). Let g be the C^'" extension of x^g-^. to M. 
According to [Le], p has the following properties. 

Lemma 2.1. p E C^'"(M) n C^'"(A'/) and has the following properties 

(a) p = p{x, y) = e^y^x + 0{x^), where v G C^'°'{dM) satisfies go = e^''^''y^x^g+{x)\j.g-j^; 

(b) p^g+ e C^'°'{M]y'^T*M) has a C^'" extension to M, such that p'^g+\rpgjj = go; 

(c) |djO|| = \ in a neighborhood of dM . 

Proof. For any q G dM, choose a local coordinate ([0, e) x U;x,y) = ([0, e) x U]y^) around q and 
write p = e^^^'^^x. Then u satisfies u{0,y) = v{y) and \d{e"x)\l2^- = 1 in [/. The constant length 
condition gives 

1 — \dx\^ 
2lf^d,u + x\du\l = 

^ X 

This equation is a non-characteristic first order PDF for u, which can be solved in a neighborhood 
of dM by Hamilton- Jacobi theory. Let F : T*M — > R denote the function 

l-\dx{p)g 

(8) F{p,0^2{dx{p),0j + ^{pM\a 

for p e M, ^ e T*M. Solving ([5]) is equivalent to finding a function u on M such that F(p, du) = 0. 

Since 'g is C'^'", the first two terms in F are C^'" functions on T*M. Letting b{p) ~ —7^^-^, 
then h e ci'"(If) n C2^"(M) and e C^^°'{dM). Thus the Hamilton vector field Xp is in 

C"''"(M) n Ci'"(M) and its flow Lp:Rx T*JT — > T*l4 is a C°'" map whenever it is defined. Let 

n n 

Lo{y) ^ du\gjY = dxu{0, y)dx + ^ dy^v{y)dy' = vo{y)dx + ^ v^{y)dy\ 

i=l i=l 

where Vi{y) = div{y). By solving ([5]) at the boundary, vo is uniquely determined and vo G C^'"{dM). 
Let ^ be the flow-out by Xp from the set {ip,Ld) & T*JI : p = (0,y) e dM}. Since F{p,uj) = 
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for p G dM, the image oi ip is a. C°'" lagrangian submanifold with the boundary of T*M contained 
in F~^; moreover, it is C"'^'" in the interior. Since it is transverse to the fibers of T*M along dM, 
so near dM, the flow out is the image of a closed 1-form uj G C°^°'(M;T*M) n C'^'°'{M;T*M). So 
locally, Ld = du for some function u e C^'"(Af ) n C^'"(Af). With u(0, y) = v{y) on the boundary, u 
is uniquely determined and clearly p — e"*^^'^^a; satisfies (a), [h), (c). □ 

Proposition 2.1. There exists a C^'" diffeomorphism $ : M > M satisfying ~ Id, with a 

C^'" restriction ^\m ■ M — > M , such that in a neighborhood of the boundary, 

(9) g+ = ^*g+=x-\dx^ +g{x)) 

where g{x) is smooth family of metrics on dM parametrized by x with Taylor expansion even up to 
order n, i.e. 

(10) g{x) = .90 + x^g2 + ■■■ + x"-^gn-i + x"g„ + 0{x"+^). 

Here g^ is the smooth representative of x'^g^\rpgjj as above and gi G C°°{dM;S^T*dM) for i = 
2, 4, n — 1 and n. 

Proof. According to [CDLS], there exist a C^'" collar diffeomorphism ^ from [0,e) x dM to a 
neighborhood of dAI, denote as W such that '^~^{x,y) = {p{x,y),y) and 'i*g+ = p~'^{dp^ + G{p)), 
where G e C°^{[0,e), Met°^{dM)) and G(0) = 50. Then ^ is C^-" away from the boundary since 
p G C^'"(M)nC^'"(A/). On the other hand, using local coordinates {x, y), we can identify [0, e) x dM 
with a neighborhood of dM denoted as W. Hence 5* can be viewed as a C^'" diffeomorphism from 
W to W, such that 'i'lgjj = Id. Let ^ be a C°° function on M such that fj,{x) = 1 is a; < -j, p.{x) — 

if a; > i and -8 < fi' < 0. For < e' < e and e small enough, -.W D [0, e'] x dM — >C W is 
well defined by 

— X 

Then for < a; < e' and e' small, dooipe') > 0. Hence : [0, e'] x dM — > [0, e'] x is a ° 
diffeomorphism, which is C^'" on (0, e'] x dM and satisfies 

[0,^]xdM 'l^,e']xdM 

Therefor ^'e' can be extend to a C^'" diffeomorphism on M. Let $ = It is easy to check that 

such $ has all the above properties. □ 

Now let's make clear the assumptions and notations which are fixed from Section [2] to Section [H 
Assume: 

(Al) For s <E {—9,6), 9 > 0, g+{s) is a family of C^'" conformally compact Poincare-Einstein 
metrics on M with a family of conformal infinities, i.e, 3 go{s) G C^{{—9, 9), Met°° {dM)) 
such that go{s) G [2:^(7+ (s)|t9m] for each s; 

(A2) g+{0) is G°° conformally compact. 

Under this two assumptions, using Proposition 12.11 there exists a G^ family of C^'" diffeomor- 
phisms over M, which are C^'" on M, such that g+{s) = ^*sg+{s) is also a G"^ family of 
G^'°' conformally compact Poincare-Einstein metrics, which are smooth near the bondary, i.e. for 
Q <x < e, 

~g+{s) = x''^{dx^ +go{s)+x^g2{x) + • • • + x"-\9„_i(s) + x"(7„(a;) + 0(x"+i)). 
In particular, $0 is smooth. Hence g+(0) is smooth. For simplicity, let 

g = .g+(0) = x-'^{dx^ +ga + x^g2 + ■■■ + x^-^gn-i + x^-g,, + 0{x"'+^)), 

h = ^^^[=0 - x-'{g', + x'g', + ■■■ + + x^g'^ + 0{x-+')). 

where g'^ — ^^%^|s=o for i = 0,2, n. Then the linearized Einstein equation 



(11) -Agh- S*gPgh + Rh + nh^O 
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is obtained by calculating j^\B=o{Ric{g{s)) + ng{s)) — 0. See [GL] for details. Here Ag, Sg and 
Sg are the standard connection Laplacian, divergence operator and its adjoint with respect to the 
metric g and 

Pgh — Sgh + —dtTgh 

is the Bianchi operator; R denotes the the order operator induced by the curvature of g. For 
h e C^{M;S^T*M) and u G C^{M;T*M), in any local coordinate patch {W C M; {y^joK^Kn), 
these operators can be written as 



The linearized Dirichlet-to-Neumann map, if well-defined at go, is the map (?q i— > 5^. We want to 
obtain this map by studying the properties of the solution to (fTTj) . The first well-known difficulty 
is that (jlip is not elliptic, but only transverse elliptic, since Einstein equation is invariant under 
diffeomorphisms. To use elliptic theory, we first need to choose a gauge to make (fTTj) elliptic. The 
condition we impose is Pgh = 0. We wish to find a family of diffeomorphisms G Diff^(M), 
such that h = ^|s=o*&s5+(s) satisfies Pgh = 0. Then h will satisfy the elliptic equation 

(12) Agh + 2Rh + 2nh^0. 

Since we are dealing with the linearized gauge condition, instead of finding exact $3 for s € {—0, 9), 
we only need to look for a generating vector field X , with dual 1-form lu, such that 

h — Lxg + h = S'gUj + h 

satisfies Pgh = 0. Then uj should satisfy the equation 

2(3gS*gUJ + 2(3gh = ■i=> AgUJ -RLJ + 2(3gh = 0. 

Since g is a Poincare-Einstein metric implies that Ruj — —nui, the linearized gauge fixing problem 
reduces to finding a solution to 

(13) AgW + nw + 2l3gh = 0. 

Set Jg = Ag + n e Diff2(M; T*M) and Lg = Ag + 2R + 2n e mS^{M; S^T*M). They are both 
uniformly degenerate elliptic operators. Next, we will carefully study these two operators, as well 
as the solutions to and p^ . The main technique is to use the 0-calculus due to Mazzeo and 
Melrose, which will be introduced in the following section. 

3. 0-Calculus 

In this section, we review the construction of a parametrix for an elliptic differential operator on 
an asymptocially hyperbolic space due to Mazzeo and Melrose ([MM]). As in Section [2l {M,g) is 
a Riemannian manifold of dimension n + 1 with smooth boundary defining function x. Near the 
boundary, 

g = x'^idx^ +go + x^g2 + ■■■ + a;""^(/„-i + a;"ff„ + 0{x''+^)). 

Hence all sectional curvatures have limit — Ma^|^2g|g-p' = —1, i.e., (M, g) is asymptotically hyperbolic. 
Suppose y is a smooth vector bundle over M with fibre dimension N. 



\A h\ = -h " 
[5gh]y = -h^^^", 

— Ra^u^h ^ -f- 2{R ^ho^y -j- R jyhfj,a), 

[Aguj]^ = -u;^./, 
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3.1. Vector Fields and Differential Operators. For each q G dM , lot {U: y) be a local coordinate 
patch in dM around q. Then ([0, e) x U; x, y) is a local coordinate patch in M around q. 

Let be the Lie algebra consisting of vector fields which are locally sums of C°° function 
multiples of the vector field xdx and xdy^ . can be viewed as the set of C°° sections of a natural 
vector bundle over M, denoted by °TM. Obviously, {xdx, xdy^ , xdy^ } gives a local basis of ^TM. 
The dual bundle, naturally denoted by °T*M, has a local basis ^}. Moreover, ^TM 

and ^T*M can be viewed as vector bundles over M, which are isomorphic to TM and T*M over 
M . However, the isomorphisms arc not natural at dM . Similarly, we can define "Pb bo the space of 
C°° vector fields generated by {xdx,dyi, dy^}. Notice that the Lie algebras and satisfy 

[%,%]^%, [%,%]<^%, 

where [•, •] is the Lie bracket. 

For any constant m € R, denote by S"'{"T*M) the space of u G C°°("T*M) such that in any 
local coordinates {W C M;p) away from the boundary, i.e. W (1 dM = 0, denoting by ^T*M\w = 
Wp X M^"'"^ canonically, then for any (n + l)-multiple indices a, (3, 

\dp'lu{p,o\<c^0{i+\(.\r-\^^ 

for some constant Cap; and in any local coordinates ([0, e) x U;x,y);x,y) near the boundary with 
U C dM, denoting by °T*M|[o,e)x!7 = [0,e) xU x 1R"+^ canonically, then for any n-multiple indices 

a, (3 and k, I G No, 

|5,^5,"4a,^u(x,y,T,r7)| < + \r\ + 

for some constant CajSkl- Similarly, we can define S'^(^T*M;End{y)). 

For m G No, let D\S^{M-,V) be the space of differential operators of order m, which can be 
written as a finite sum of at most m-fold products of vector fields in Xo with values in End{V). 
Near the boundary, in a coordinate patch ([0, c) x U;x, y), choose a trivialization of V, such that 
V[o,e)xu = ([0, e) X 17) X C^. Then any P e Diff^(M; V) is of the form 

(14) Pij= Cij,ma{x,y){xdxr{xdyr. 

\a\+m<k 

The 0-princinpal symbol of P is defined by 

fc+|a|=m 

The coordinate invariant property shows that ^am{P) is globally well defined, i.e. 

°(Tm{P) e S'^CT*M;End{V))/S"'-\°T*M;End{V)). 

Using the symbol language, we can define the 'small' 0-pseudo- differential operators ^^{M;V) 
naturally as the usual smooth pseudo differential operators. The symbol map induces a short exact 
sequence 

(15) — > *^-i(M; V) *^(M; V) ^ S'^{°T*M; End{V))/S'^-\°T*M; End{V)) — > 0. 

3.2. Stretched Product Space. To investigate the properties of Schwartz kernel for operators in 
^'™(M), we recall the stretched product M Xq M obtained by blowing up the boundary diagonal 
diag{dM) in M x M. Let (3 be the blow down map, i.e. 

p-.MxoM — >MxM. 

If we denote by {x',y') the corresponding local coordinate on the second copy of M, then in the 
neighborhood around {q,q) in diag{dM), we can introduce the polar coordinates 

R=\x'^ +x''^ + \y-y'f\^, r = \y-y'\, s=^, t=—, 

x' X 

y-y' I y-y' x , x' y-y' 

X X K K r 

which give the local coordinates in M Xq M around (3~^{q, q). M Xo M is a compact manifold with 
corners up to dimension 3. It has three boundary hypersurfaces, top face T, bottom face B and the 
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front face F. Then R (rcsp. p, p') is the defining function for F (resp. T. B). The front face at q is 
denoted by Fq = (3~^{q, q), which is a smooth compact manifold with corners of codimension 2 and 
has two boundary hypersurfaces -Fq fl T and FgdB with corresponding defining functions p and p'. 
Moreover, the interior of Fq is isomorphism to hyperbolic space IHI"+^, or the inward half tangent 
space T^M, carrying on the nature semi-product group structure as M_|_ k R". Denote by diagQ{M) 
the closure of I3~'^ {diag{M)) in M Xq M. Clearly, 

diagoiM) f) F = {s = l,z = 0,R= 0}. 

Let [M xqM]^ be the % stretched product doubled across the front face. It has the structure of a 
smooth manifold with corners up to codimension 2. Furthermore, the doubled diagonal [diagQ{M)]'^ 
lies in the interior. Note that any vector bundle over M Xq M can be extended to [M Xq M]^. 

Similarly, dM XqM is the blow-up of dM x M along diag{dM). It can be embedded in M Xq M. 
Actually, the identification of dM Xq M with {p = 0} C M Xq M provides the embedding map. 
Denote by (3 the blown down map 

0:dMxoM — >dM xM. 

The front face of dM Xq M at any q e dJd is just FqHT. 

Finally, denote by dM Xq dM the blow up of dM x dM along diag{dM), which is naturally 
embedded into dM Xq M by identification with {p' = 0} C dM Xq M, with the blow down map 

0:dMxodM — ^ dM x dM. 

Note that dM Xq dM is a smooth manifold with boundary containing the lift of diag{dM) via (3, 
which can be identified with the spherical normal bundle of diag{dM), i.e. SN{diag{dM)). 

1 

3.3. Half Densities. Let T§ (M) be the line bundle of singular half densities on M, trivialized by 

v := \dvolg\^ , and T^i^dM) be the bundle of half density on dM, trivialized by uq := \dvolgg\^ . 

Naturally we can define T§(M x M) (resp. T^{dM x M) and T^dM x dM)) the half density 
bundle on M x M) (resp. dM x M and dM x dM) trivialized hy v x v (resp. pq x v and x vo)- 
Hence finally we can define the half density bundles over the blown up spaces via the blow down 
map, i.e 

r| (M xo M) = /3*r| (M X M), 
r| {dM xo M) = ^*r| {dM XqM), 
r|(9M xo dM) = P*T^{dM x dM). 

Note that T§ {M Xq M) is trivial on Fq. For example, in coordinates {x', y', s, z), T§ {M Xq M)\f^ 
is spanned by | ^?n+^g1f+i |^; in coordinates {x,y,t,z'), it is spanned by | ^n+i ^1+^ 1 1 ^ ; in coordinates 
{p, p',r,uj, y) , it is spanned by | ^ff ^[f^l'^f^i | ^ . Similarly F^ {dM x o M) is trivial on Fq n T. Hence, 

r|(M Xo M) extends to half density bimdle F|([M Xo M]^) over [M Xq M]^. 
With such half density bundles, we can define the Sobolev spaces by 

Lg(M;r|(M)) = {/eC-~(M;r|(M)): / |/|^ < oo} 

which is a Hilbert space with inner product 

(/1)/2)l2(M) = / fih 

Jm 

for any /i,/2 S i^(M;F|(M)). Moreover, for m e N, 

H^{M;rl{M)) = {/ e LI{M;tI{M)) : f e i^(M; r| (M)), < k < m}, 
H^^{M-Tl{M) = {H^{M-Vl(M))'. 
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3.4. Boundary Cornormal Distributions. Let X bo any compact manifold with k boundary 
hypersurfaces diX for 1 < i < fc, with corresponding defining function pi. Let if G No be a fixed 
constant. Then for any ci, ...,Cfe e C, we can define a functions space - associated to K 
with certain conormal sigularities at the boundary surfaces, i.e. 

N 

^^^^ _^c„...,c.(^) ^ 1^ g : r^u e • • • pl^{Y^{\ogp,r^^ ■ ■ ■ (logpfe)'^-L-(X)), 

for some e No and for 1 < i < TV, 1 < j < fc, € No ("1 [0, K]} 

And if some = oo, define ^ci,...oo,...,c,(j^) ^ fli^i ■'^''''■■■''■■■'''H^)- 

Moreover, we can define ■^plg {X) to be the natural subspace of polyhomogeneous elements 
in by replacing L°°{X) by C°°(X) in the definition, i.e. 

N 

(17) Kkr"^^) = e ^'(X) : Yiu G • ■■ptiY.i^ogp.f'^ ■ ■ ■ (logpO'""C^°^(X)), 

for some iV G No and for 1 < i < AT, 1 < j < k, dij G No n [0, K]} 

These definitions can be easily extended to boundary conormal distributional sections of a vector 

bundle over X. 

3.5. O-Pseudodifferential Operators. Let Jf^"^ (M XoM; diago{M); End{V)(g)r^) be the space of 

j_ 

distributional sections of £n(i(F)(g)ro (MxqM) onMxoM, which can be extended to [MxqM]'^ as 

distributional sections of End{V)®rQ ([MxqM]^), conormal along the doubled diagonal [diago{M)]^ 
of order m and vanishing to infinite order on the doubled boundary faces T and B. Then 

P G *^(M; y O r|) <^ k(P) G ^™(M Xq M; diago{M)- End{V) O r|), 

where k(P) denotes the Schwartz kernel of P. 

Suppose V\qj;j is decomposed into ® ■ ■ ■ where is a vector bundle over dM for 

1 < i < r < N . Given any constant vector a = (ai, a^)* G C. Then for some constant K G No, 
we can define 

r if 

(M; y ® r| ) = {i; G ^' (M; F r| ) n C°° (M; y ® r| ) : = ^ ^ (log a;)^' u^,- , 

(18) i=i j=o 

Uij- G C°°(M; V ® rj), Uy lgj^j e C°°(9M; T/* ® r}), for 1 < i < r and < j < K} 

Given constant matrices Et, Eb G M (r, C). Let [ET\i, denote the the i-th row of Et and [-Eb].^ 

denote the j-th column of Eb- Then we can define s^phg^'^iW Xo Idf-.Endiy) ® r|) to be a 

subspace of ^'([M Xq M]^; i;nd(T/) ® r|) such that any element u e ^/^^^'-^^ ([M Xq M]^; £;nd(y) ® 

Fq ) has the form [u]»j ~ ubp'^^^^'^ near boundary face i? and [iijj, ~ pl^^l^'ur near boundary face 
T, where with some cut off function x satisfying x(p) = 1 if |p| < e and x(p) = if \p\ > 2e for 

some c > 0, x{p')ub, x{p)ut & £/°'°{[M Xo M]^); End{V) (g) T§). Note, for any constant c, and any 
defining function x, we always denote 

[Et + c],j = [EtI,, + c, = xl^^l- . 

Now we can define the space of distributions cornormal to boundaries T and B and smooth up 
to the front face F and the corresponding space of operators: 

(19) jro-°°'^-'^-(MxoM;i;nd(F)®r|) =^^-'^«([MxoM]^;^^^^ 

*o°°'^^'^^(M;y (g)r|) = {A € C-°°{M;End{V) <»rl) : 

k{A) g j(f-^^i^T,EB Xo M; End{V) (g) r|). 
So we have the 'large' class of 0-pseudo- differential operators defined by 
(21) *™'-^^"^^(M;F(8)r|) = 4'™(M;For|) +*-°°'^^'^«(M;F(g)r|) 
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with corresponding kernel space 

j^n^,ET,EB^j^ XfjM;diagoiM); End{V) (^rl) = ^^^^"(M XqM; diago{M); EndiV) ® r|) 

^ .^ocE^.Es (7g X M; End{V) ® T J ) . 

Note that the decomposition of V, which has only been considered at the boundary so far, usually 
can be extended to a neighborhood of DM up to 0{x^) for some certain constant A e No U {00} in 
practice. Suppose this is the case. Consider u as a section of V. Near the boundary, if writing 

r K 

u = ^^a:'^'(logx)^M,„with G C°°(M;F0rJ)), u,,\g^ £ (^M; ® r| ) , 

i=l j=0 

then the expression makes sense if and if only 

\^a,-^aj\<A, yi<i,j<r. 

Usually, A is determined by V, M and the metric. Hence ^f™'^^'^^ (M; T^^Fq) makes sense if and 
only if Et, Eb satisfy the following condition: 

Condition 1. 

\min i<j<r{^[ET]ij} - min i<j<r{^[ET]k]}\ < A, y I < i,k < r, 
\min i<j<r{^[EB]ji} - min i<j<n{^[EB]jk}\ < A, V 1 < z, fc < r. 

According to [Ma], we have the mapping properties for such large class of 0-pseudodifferential 
operators. 

Proposition 3.1. For each A e ^q''^'^'^'' {M;V (g) r|) and constants m, m' e Nq and b, b' e M, 

A : x''H^'+"''{M; F (g) tJ) — > x'''hS''{M; F (g) r|) 
is a bounded operator provided that 



min {[Erhj} -b' > min {{Eb^j} + b>-, 



min {[ET]rj} + min {[-Bsly} > n, b> b' . 
Proposition 3.2. For any A e qi'^'^^T^En ^j^j. ^/ ^ r|), let e-j = [Et],j for l<j<r. Then 

r 

A : C'°°(M;F(g)r|) — > ^ (M; F (g) r|) 

is continuous. 

3.6. Normal Operator. Suppose A E ^fj^'^^'^^ (M; F (g Fq ). The normal operator of A at g is 
defined by 

N{A)f^ [ n{A)\FJ{-,y--z)-dz■d^l, V/ e (H"+i, C^), 
J s s s 

where \ . Note that F|f, = 7r*(V^|5) is trivial, where tt is the projection map T+Af {q}. 

In particular, if P e DiHo"(Af ; V ®Vq) is of the form (fT4|) . then its normal operator is 

\oL\+m<k 

The normal operator characterizes the behavior of the original operator near the boundary, which 
can be easily understood via an equivalent definition 

Nq{P)u = lim i?;$*P($"i)*i?lu, 
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where $ is a local diffeomorphism carrying a neighborhood of G H"+'^ to a neighborhood of g G dM 
and Rt denotes the scaling. According to [MM], we have a short exact sequence 

^™,ij,,nT,B.,nB ^p^, ^ dtagoiM); M {N, C) ® Tj) ^ 0. 

Here R is the boundary defining function of Fg in U Xq IT and i?*™"^^'-^^ (M; Snd(F) (g) T^) 
consists of operators with kernels in (M Xq M; diago{M); End{V) Fq ). 

3.7. Indicial Operator. Also for A e ^'""^ (M; F r|) we can define the indicial operator 
by 

IiA){s)f = x-%A{x'f))\g-j^, V/ G C°°(9M;F® r|). 
Similarly, if P G Diffo^(A/; V <^T§) with components as in (fT4|l . then we can write 

= E C»,,™o(0,y)s"\ V,s G C. 

The indicial operator characterizes the leading order of the operator near the boundary. According 
to [MM], there is a short exact sequence 

(24) 

— > i?°°*o°°'-^^+^^-^^(Af;£;nd(F) ® rj) — > R'^-^-°°''^^'^'' {M; End{V) (^Tl) 
R^.je^~°°'^^'^''(M XoM;EndiV) ® r|)/i?°°jro"°°'^^+^'^^ (M XoM;EndiV)®rl) — > 0, 
where R is the boundary defining function of the front face F in M Xq M. 

3.8. Parametrix. The construction of a parametrix for an elliptic differential operator of 0-type 
is based on the three short exact sequences (|15p . (1^51) and (|24l) generated by symbol map, normal 
operator and indicial operator respectively. Here we only do this for a small class of elliptic differential 
operators, which includes Lg and Jg introduced in Section [2l First, notice that 

Lemma 3.1. If A £ (M; «) r| ), B G R°°%°°-'^^''^'' {M; End{V) (g> tI) and 

Pg Diffo(M;i;nd(T/)(8)r|), then 

N{PA) = N{P)N{A), 
I{PB) = I(P)I(B). 

Theorem 3.1. Let V he a vector bundle over M with fibre dimension N and a decomposition 
^laM ~ ® ■ ■ ■ ® y at the boundary, for some 1 < r < N . Suppose that the decomposition 
can be extended to a neiborhood of dM up to 0{x^) for some A G N U {oo} and suppose that 

P G Diffo(M;y (K) r|) is elliptic, i.e., °a2{P) G S'^[^T*M;End{V) ®T^) is invertible. Assume P 
also has the following properties: 

(1) The indicial operator I{P) is diagonal. For 1 < k < r, I{P){s)\k ~ fk{s), where fk{s) ~ 
— + ns + Ck for some Ck > 0. Let Sk < s'' be the two indicial roots of fk{s) — 0, then 
\si ~ Sj\ < A for I < i, j < r. Moreover, 

r 

'-^^{s,s)^f]{s,,s'^)^il}. 

(2) Vg G dM, the normal operator Nq{P) is invertible, i.e. 

NgiP)-' : C^iFg-.MiN, C) ® Tj) ^Z'^'^iFg; M{N, C) ® r|) 
is bounded, where Et and Eb are r x r matrices satisfying Condition Q] and 

[Erjij — + I — 6ij, Eb — E^. 
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Then there exists a right parametrix Q G ^'q '"'^^^^^^ (^]\,f-^ End{V) (8" Fq ) such that 
PQ-Id = Rie R°^^^^'°°''^''{M-End{V) ® T^). 

Proof. Use the symbol calculus and (jlSp . it is easy to find a first approximation to the parametrix, 
Qi e *(7™(M; End{V) (g) r| ) such that 

PQi -Id = Eie *o °°(M; End{V) ® rj). 

Notice that the Schwartz kernel k{Ei) vanishes to infinite order at the boundary faces T and B. 
Hence 

k{Ei) r-^^R^Kj, e C'°°(B"+i;M(iV,C)®r|) e L^(B"+i;M(iV,C)®r|). 

j 

Since N{P) is invertible, by there exists Q2 G \^^°°-^t,Eb g^^^j^ ^-j-^g^^- 

PQ2-Ei^E2 e R°°^^°°'''^'''"{M;End{V)®Tl). 

Notice that P actually kills the leading term of Q2 at boundary T, by a modification as in [MM], 
we can choose Q2 e xi,-'^-^^T+id..Es jjg^^g 

PQ2-Ei^E2 e R°^%°°^^^+''^'^^ (M; End{V) ®tI). 
So k{E2) e ^Et+m.Eb^ At last, since I{P){Et + Id + c) is invertible for aU c e No, we can find 
Q3 e R'=^%°°-^^+''^-^''{M;End{V)(g)T^), such that 

PQ3 - i?2 - ^3 e ® r|). 

Then Q = Qi +Q2 + ^3- □ 

4. Gauge Fixing 

In this section, let's study the operator Jg introduced in Section [5] carefully. Recall near the 
boundary 

g = x-^{dx^ +go + x^g2 + ■■■ + x"-^gn^i + x^-g,, + 0(a;"+^)), 
and Jg = Ag + n acting on sections of T* M . More naturally, we can view Jg as an operator acting 
on sections of °T*M. 

4.1. Indicial Operator. Clearly, ^T*M\qjj has a decomposition with repect to metric i.e. 

Such decomposition can be extended to a neighborhood of the boundary by identifying it with 
[0, e) X dM in the geodesic normal way, i.e. 

where the ' _L' is defined by the induced metric on ^T*M from g. Hence in this case, we can choose 
A = 00 in Condition [1] The indicial operator corresponding to this decomposition is 

2n 
(n+1) 

So the two indicial functions corresponding to this decomposition are 

/i(s) = — + ns + 2n and /2(s) = —x'^ + ns + n + 1 
with corresponding indicial roots 



I{P){s) = -s^ +ns 



V s e 



n — \/n'^ + 8n -1 n + + 8n , 
si = , s — and S2 = — 1, s — n + 1. 



Hence (s, s) = (—1, n + 1). 



12 



FANG WANG 



4.2. Normal Operator. To compute the normal operator, consider the half plane model of hyper- 
bolic space 

H"+i=M+xM" with metric ^ = ^^!±^^ dvolr = ^. 

For any q e dM, °T*M has a basis {^,^,...,^}, which corresponds to the basis of "T*H"+^ 

Set e° = ^ and = ^ for 1 < i < n. Then Ot*H"+i = e"+i x (C C") corresponds to the 
decomposition of ^T*Mgj^. The normal operator of Jg can be expressed as 



i=l 



2n 
(n + l)Idn 



+ 2x 



-dyl 
dyl 



L 5,. 



On the other hand Ng{Jg) = = + n acting on °T*H"+^, so it is self-adjoint. 
Lemma 4.1. J'or any a e ( "~^^'+^" , "+^^'+^" ), 

is an isomorphism. 

Proof. For each u e F^(H"+i; C"+i), JVq( Jg)(a;"n) = a(n - a)x''u + x''Ng{Jg)u. Hence, 

{Ng{Jg){x''u),x"u)^aL^ = {a{n- a)u + Ng{Jg)u,u)L2 > ||Vn||^2 + {-a^ + na + n)\u\l2 > 0, 
and equality holds if and only if u = 0. □ 

For any u e Lg(H"+i, C"+i), let u{x,£,) = .^y^^(u)(x, ^) be the Fourier Transform with re- 
spect to the variable y. Then Ng{Jg)u{x,y) = is reduced to the ordinary differential equation 
Nq{Jg)u{x,^) = 0, where 



2n 
{n + l)Idn 



+ 2ix 



-e 
e 



C 











Letting x = x\^\ for ^ ^ and ^ = ^/|^| G S" ^, then the equation Nq{Jg)v{x\^\,£,) = is equivalent 
to Ng{Jg)v{x,^) = 0, where 



2n 
{n + l)Idn 



+ 2ix 



-i^ 





Ng{Jg) = -{xdsf + nxds + x^ + 

L 

Let r|(M+) be the half density bundle over R+ trivialized by l^^l^. 
Lemma 4.2. For any ^ € 5""^, the map 

is continuous and positive. 







i^(K+;C"+l®r| 
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Proof. For any w{x) = J2^=o ^^^^ ^ -f^o (^J ; C"+^ ® 

71. 

{Nq{Jg)w,w)L2 = \\xdsw\\l2 + \\xw\\L2+272\\wo\\L2 + {n+l)J2\\wi\\Ll 

1=1 

« n 

Jr+ 

> \\^dxw\\li + + MwoWl^ + {n + 1) ^ Ik.lLg - \\w\\l2 - \\xw\\l2 

i=l 

> \\xdsw\\l2+n\\w\\L2 > 

And equality holds if and only if w = 0. □ 



Lemma 4.3. The map 

is continuous and positive 



N,{Jg) : H^{Rt; {0-'L\R'l, C"+i) ® r| (R+)) — > ^^(Kn. ^ r|(M+)) 



Proo/. For any w{x,^) e H^{R+; {O^^L^{R'';C''+'^) (E) r§{R+)), by changing variables and last 
Lemma, we have 

{Nq{Jg)w,w)L2L^ > \\xd^w\\l2^2 +n\\w\\L2L^ > 
and {Nq{Jg)w,'w)i^2j^2 — holds if and only if w = 0. □ 

According to [CL], the linear ODE Nq{Jg)u = has a regular singular point at a: = and irregular 
singular point at ai = c» and no other singular points. So a formal solution at x = is a true solution 
and for each formal solution at a; = oo, there is an actual solution which has asymptotic expansions 
at a; = oo corresponding to it. Hence, to find all 2N linear independent solutions of Nq{Jg)u = 0, we 
only need to find out all the linear independent formal solutions at a; = 0. Since such solutions can 
not blow up at any point between and oo, we can extend them to the whole half real line (0, oo). 
By using the indicial functions fi and indicial roots Si, defined in Subsection 4.1 for i = 1, 2, the 
construction of formal solutions at a; = is simple as follows. 

(1) Start with a'' e°. Since fi{s^ + k) ^ Q for all /c g N and i = 1,2, so assume w° = a* e" + 

^^Wk is a solution to Nq{Jg)u = we can solve each function Wk by induction. 

(2) Start with a^'^e'^. Since fi{si + fc) 7^ for all fc G N and i = 1, 2, as in (1), there is a solution 

(3) For each 1 < i < n, start with a* e*. Since /i(s^ + A:) ^ for all fc e N and i = 1, 2, as in 
(1), we have a solution of the form a* + YlV=i ^^'wk- 

(4) For each 1 < « < n, start with x'^^e*. Note in this case, /i(s2 + fc) 7^ for all fc S N and 
/2(s2 + fc) 7^ only if fc 7^ n + 2. Write v^ = a^^e* + J2T=i x'^^'^Wk + x"+^ log aw, we still 
can solve all Wk and w by induction to makt Vi a solution. 

On the other hand, the formal solution at a: = 00 has leading term either x^e^^ or x^e^. By 
Lemma |4.2[ as a ^ 00, for < i < n, the behavior of solutions constructed above is 



Since Nq{Jg) is self djoint with respect to a ^'^'^^^dx, the standard method in [CL] yields the 
Green function for Nq{Jg), i.e. 

Gi(a, x', i) ^ U{x, i)V*{x\ i)H{x' - a) + V{x, i)U*{x\ i)H{x - a'), 

where J7(a,^) and V{x,£,) are (n + 1) x (n + 1) matrices whose columns are linear combination of 
and Vi. 

Nq{Jg)Gi{x,x' ,i) ^ x^'+Hix - x'). 

Moreover, 

U'{x,i)V*{x',i) - V'{x,i)U*{x'A) = -ix"-i/ci. 
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This implies that 

U{x,^) = (aot;°,aiw\ ...,a„i;"), V{x,^) = (bovo^bivi, ...,6„i'„), 
for some ^ 0, bi ^ ioi < i < n. 

Lemma 4.4. The Green function Gi{x,x' induces a bounded integral operator on L'^{M.'^ (M.^)) . 
Lemma 4.5. There exists a unique kernel G2{x, x' , £,) which is bounded on Lq(M+; L^(]R^; C"^"'^) (g) 
rJ(R+)) satisfying 

Nq{Jg)G2{x, x', = x"+^6{x - x') (g) Id. 

Proof. Let 

G2(x,a;',0 = |er"Gi(a:|e|,a;'|e|,|) 
Since the 1-dimensional 6 function is of homogeneity — 1 , it is easy to check that 

Nq{Jg)G2{x, x', i) = - x') (g) Id. 

Note for any a > 0, the map 

fi^) — ' fa{x) = a'^fia-^x) 
is an isometry on ig(R+; r| (R+)). Therefore, for any u{x,Cj e Lg(M+; ^^(M"; C"+i) (g) Tj (M+)), 



G2{x,x,£,)u{x',£,)) 



dx' 



= / / 



G,{x,x\m\^<\^\-'^',o) 



dx 
dx' 



dx 

rn+l 



< cm\^um-'x,o\\h^4u\\] 



Integrating the above inequaUty over gives the Lg(H"+-'^; C"'^"'^ Fq )-boundness of G2(x, x' ,£,). 
The uniqueness follows directly from Lemma [4.31 i-e- there is no Lq null space. □ 

Proposition 4.1. The operator Nq{Jg) = Jr : i7^(H"+i; C"+i r|) — > L2(H"+1;C"+i (g) r|) is 
invertible with Green function 

1 



(25) 



G{x,x',y,y') = 



(27r)" 



|^r"Gi(x|e|,x'|^UV«(^-^')de 



Moreover, G £ ^-^.-^t.-Eb (h»+1; M(ri + 1, C) r= ), where 



2 



n + 2 



2 2 

n+2 n+ 1 



Proof. The green function constructed as ([25]) satisfies 

iV,( J<,)G(x, a;', y') ^ - x')<5(y - y'). 

The Lg(H"+i; °r*H"+i (g) r|)-boundness of G follows directly from Lemmag^l The behavior of G 
at the boundary T and B is obtained by computing the Fourier Transform (I25p based on asymptotic 
analysis of G2{x, x' ,£,). We can follow the proof given by Mazzeo in [Ma] step by step. □ 

Notice that the front face can be thought of as the ball B"^^ blown up a point on the boundary. 
We can change all the discussion above for EI"+^ to the ball model Let ei = (s^, s^ + 1)* and 

62 = (s^ + 1, s^)* corresponding to the i-th column of Et for i = 1, 2. 



Corollary 4.1. The map J-^ : C'°°(s"'^\ C"+i (g Lq^ ) — > ELi •</Jg(^"^'; (g Lq^ ) is con- 
tinuous. 

Identifying Fq with _B"+^ blown up p e where p corresponds to the original defining 

function of _B"+^ and p' is corresponds to the new defining function generated by the blow-up. A 
similar analysis as in the proof of Proposition 6.19 in [MM] shows that 



-n+l 
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Corollary 4.2. The map Ng{Jg)-^ : ^""(T^; M(n+ 1, C) (g) rj) — > (T^; Af (n+ 1, C) (g)r^) 

is continuous. 

Notice that the foundamental solution U{x, ^) does not involve any logrithm terms. Hence in this 
case, we can take the constant K — in the definition of pT|. i.e. all the spaces ■B^pi^g{M) used in 
this section can actually be replaced by C°°(M). 

4.3. Parametrlx. 

Proposition 4.2. The operator Jg : H^{M-"T*M (g) r|) — > Ll{M;°T*M (g) r|) is self adjotnt, 
positive and thus invertible. 

Proof. Clearly Jg ^ V*V + n is self-adjoint. For any uj G H^(M;°T*M (g) r|), 

{Jg^.^)x-Ll = l|Vw|li2 +n\\uj\\Li > 0, 
with equality holds if and only if a; = 0. So Jg is positive and thus injective. Note Jg satisfies all 
the assumption in Theorem 13. li so there exists Q G \I>^2,-Et,£b ^^j.^ '^T*M ®Tq), such that 

JgQ-Id = Ri e i?°°*o°°^°°-^^(M;°T*M(gr|), 
Q*Jg-Id = Rl e i?°°^'o°°'^^'°°(M;°r*Af (grj), 

where 

Et = 2 " + ^ 

By Proposition 13. 1[ the map 

Q : Hl{M- °T*M «) r|) — ^ LI{M; °T*M (g) r|) 
is bounded. It is also clear that the two maps 

i?i : ig(A/;°r*A/(gr^) — » C'°°(Af; °T* Af (g T^) 

2 

Rl : l2(M; °T*M r|) ^ <,;,(M; "T*A// ® r|) 

1 

are bounded. Therefore, as a map from Lq(A/; ^T*M (g Fq ) to itself, i?i and R^ are both compact, 
which implies that 

Id + Ri: Ll{M; °T*M (g r|) — > LI{M; °T* M (g tJ) 

is Fredholm, so Ran{Jg) is closed in Ll{M;'^T* M) with Coker{Jg) of finite dimension. Since Jg is 
self-adjoint, Coker{Jg) — Ker{Jg) — {0}. Hence, Jg is surjective and thus invertible. □ 

Proposition 4.3. For N > 2n, Jg^ : £!/°(M;°T* M (E) T^) — > J2'i=i-^'''(M;°T*M (^T^) is 
continuous. 

Proof. Since J^"^ = Q* +RlJg'^ , the continuity of Jg'^ is directly from the continuity of the following 
map: 

2 

Q* : x^^°(M;°r*Af (gr|) — > ^ ^"^^ (M; °T*A'/ » r| ), 

1=1 

Jg^ : a;^^°(M;°r*A//g)r|) — > ig(Af ; "T* A/ ® r| ), 

2 

Rl : Ll{M,^T*M ®vl) — > ^ jz/'^' (Af; °T* AT (g r|). 

i=l 

□ 

Proposition 4.4. J^^ e ^^'^^'^'^ {M;°T*M (E)T^). 



Eh — EX. — 



n+^/n'^ +8n n+^/n'^ +8n 
2 2 

n + 2 n + 1 
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Proof. Use the notation as above, we have 

J-i =Q + J-^Rx = Q* + RIJ-^ =Q*+ RIQ + RU'^Ri = Q + Q*Ri + RU'^Ri. 

Notice that k{Ri) e i?°°=^-°°'°°'-®s(M Xq M;End{v) O r|) and k{Ri) e R°°J(r-°°'^-^'°°{M Xq 
1 

M; End{v) ® Fq ) and both vanish to infinite order at the front face. A theory of composition of 
0-pseudodifferential operators shows that 

k,{Q*Ri),k{RIJ-^Ri) g i?°°=^-°°'-^^'-^s(M XoM;End{V)®T^). 
The conclusion follows directly. □ 

4.4. Preliminary Computation. Let ^ = x^g. Then 'g\rpQjj = go- Let T\j (resp. T''^j and T'^^j 
) be the Christoffel symbols for g, (resp. 'g and go), and let Vj, (resp. Vj and Vj) be the Levi- 
Civita connection with respect to g (resp. g and go). Throughout this section, let E{z;c\, ..ci) be 
a polynomial in z without constant term and with coefficients being functions of (ci, ...,c;). Let T; 
be a differential operator of order < I on dM with coefficients being functions of the components 
of go and their tangential derivatives; Let Ti be a differential operator of order < I on dM with 
coefficients being functions of the components of go,gn and their tangential derivatives. Let T/ = 
and = if ; < 0. Recall that 

h = x-\g'o + x'g'^ + • • • + x-'g'^_, + x^g'^ + 0{x"+')). 
Lemma 4.6. There exist differential operators Ti,To such that trg^g^^ — 70.90' ^go9n — ^i^o- 
Proof. According to [FG], trg^(^s)9n{s) = and Sgg(s)gn{s) = for all s £ {—9,9), which obviously 
implies that trg„g'^^ = -trg^g'„ = Tog'o and 5g„g',^ Ti.g^. □ 

According to [FG], the one to one correspondence of normal forms for Poincare- Einstein met- 
rics and normal forms for ambient metrics implies that in local coordinates {U;y^ , ...,y'^) on the 
boundary, 

[92{s)]ij = -Pij{s) = -]^{Rij{s) - 2§^[50(s)]y), 

and for 2 < ; < 

[9^ii^)h = („-4)(n-6)".'..(n-2o (^ro(/)-P^^- " A-(4^V,^V,ir,„(,)P(s)) + Zota, 

where Rij{s) and S{s) are the components of the Ricci curvature and scalar curvature of .9o(s) and 
*V is the covariant derivative with respect to go{s). Here lots denotes quadratic and higher terms 
involving fewer derivatives of curvature. 

Lemma 4.7. By direct computation, we have 

(a) For2<l< 

92 = -7^{l^9o9o - ^goPgo9o) + 2(n-i)(n-2) (Ago^?-gogo " trg„5gJl^g'o)go + Txg'o 

921 = - 2'-l(n-4)-(n-20 {^J^d^go^'o " ^ ^*gol^3o9'o) " 2(n-2)(n-l) i^''go^'^3o9'o 

-A'g-Hrgj;jgJo)9 ~ V V ( A^; ^ frg^ ff^, - A'g-Hrg,6;jgJo)} + ^bi-i^o 

(b) IftrgJo = 0,for2<l<^ 

92 = -TiradAgo^O - ^*go^go9'o) - 2(«-l)(n-2) ft'^go'^go'^gogo)go + Tig'o 

921 = ~2'-i(n-4)---(n-2i){;;^(^Ago5o - Ag^^ Sg^Sg^go) + 2(n-2) (n- 1) ( ^go ^ *^9o '^go '^9ofl'o)fl' 

+ 2i;h)^^iKfirASg,g'o)} + T2i-i9'o 

(c) Iftrg„g'o = dgJo = 0, /or 1 < / < 

52i = d2iAg^gQ + T2i-igo 

where daj = - 2'(n-2)(n-4)...(n-20 " 
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Moreover, for 1 < I < -^^y^, each differential operator T21-1 appearing above has coefficients depend- 
ing on the components of go and their derivatives up to order 21. 

For any H G C°°((0,e) x U;S^T^M)), denote T^^^iH) = ^{da^Hp^ + dpH^p ~ d^H^p) for 
simplicity. 

Lemma 4.8. Suppose that ((0, e) x U; x,y^ , ...,y") is a local coordinate system near the boundary. 
Set X = y^ . Then the Christofel symbols of g are 

r°oo - rOo. = ^0 - r'oo - for^>l, 

= x-\~5^ + E{x^;go,...,gn-i) + ^x^'[go]'''[9nh+Oix''+')) 
for j,k> 1, 

r% = x-\[goh + E{x^-go, ...,g„_i) - i^x"[5„],, + 0(x"+i)) 
for i,j>l 

T% = f%.+i?(x2;go,...,<?„-i) + ^"A%.(go,<?n) + 0(x"+i) 
for i,j,k > 1. 

where A\j{go,gn) = [go]''*T,jt{9n) - [9nf^Tzjt{9o) = ^gof^^tignh + Vj[.9„]»i - Vt[g„]jj). 
Lemma 4.9. In local coordinates {{Q , e) x U ] x , y) , suppose Pgh — [Pgh]Qdx + '^"^^[Pgh]idy'^ . Then 

n-l 

[Pgh]o = ^""'((^ - ^>9o92i + T2i-ig'o) + x^'-'Tog'o + 0(a;"), 

1=0 

[Pghl = J2 ^'\[6aog2ih + ^^(.trgJ2l) + T2ig'o) + ^''T,g', + 0(x"+i), 
for 1 < i < n. In particular, iftig^gQ = Sg^gQ — 0, then the two formulae reduce to 

2 ^ 

[Pgh]o = J2 ^''^'(T^2i-i.9o) + ^""'^0.9^, + 0(x"), 

1=0 
n-l 

[Pgh], = ^''(^2;5o) + x^'T.g'^ + 0(x"+i). 
;=o 

Note for Q < k < n — 1, the coefficients of each differential operator Tk appearing above only depend 
on the components of go and their derivatives up to order k + I; and for k — 0,1, the coefficients 
of each differential operator Tk only depend on the compeoents of go,gn o,nd their derivatives up to 
order k. 

Proof. Note that Iiq^ — /i^o = for < /x < n. By direct computation, 

T^-l 
2 

trgoh = 2^''Ko52/ + T2i-i9'o) + ^^"(trgoff; - trg„.9^,) + 0(x"+i), 

1=0 

n-l 

[5gh]o = iZ ^''"'(~t^9o52i + T2i-ig'o) + x^'^\^tVgM + 0(x"), 

(=0 

andfor \ < i < n, 

n-l 

[Sgh], = Y ^''(^921 + T2ig'o) + ^''{[Sgognh " ['^S^Sol. 
1=0 

+ [9oy'iA%,i9o,9n)[go]sk+A',ki9o,9n)[9oh)) + Oix"+'). 
Using f3gh — Sgh + ^dti-gh and by Lemma [4.61 we can get the formulae easily. □ 
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4.5. Choice of Gauge. 

Proposition 4.5. There exist a family of C"~^ diffeomorphisms s '■ M — > M, such that 
h — ^|s=o(^s^+(s)) satisfies (3gh — 0, where s is the restriction of'^g to M. Moreover, h G 
C"'°'{AI; S'^C^T* M)) for any < a < 1, with asymptotic expansion 

h = x^^iho + xhi + ... + x^'hn + \ogxH + C'(a;"+2 logx) 

at the boundary, where Hq = gp — ^tfgod'n: ~ 2^i.9o ^ < i < n~ \ and hn — g'^ — Tng^. Moreover, 
ift'Tgag'f) = 5ga9o = 0, then ha = g'„ and /i„ = g[-^ - T„_ig^. 

Proof. As the discussion in Section [2l we can assume that X is the infinitesimal generator of the 
curve of diffeomorphisims ^'s- Then 



h^-U=o{^^G+is)) = Lxg + h. 

So we only need to show that there exists a vector field X on M such that 

Pg{Lxg + h) = 0. 

Denoting by uj the dual form of X, then Lxg = S*oj, so ui should satisfy the equaiton 

Pg{S;u;)+l3gh = 0. 

This is equivalent to the equation 

(26) Aguj - Ricgicu) + 2Pgh = 0. 

where Jg^ Ag- Ricg = Ag + ne Diffo(M; °T*M) studied in the previous subsections. To simplify 
the following computation, for any q G dM, choose normal coordinates {U; y) around p with respect 
to 5o- So at p, 



Jg = -{xdxf + nxdx - X! ^v' 



2n 

{n + l)Idn 



2x 



-dy 
dyi 



dy. 



'dy 





xD, 



for some D e Difro(M; °T*M). For oj G C°°{dM; °T*M), write lu = wq^ + ^"^^ w^^. Then 
Jgix'^Lu) = a;'=(/(Jg)(fc)u;)+a;'=+i(-^2a,u;,— + 29,u.o— +Tow) 

2—1 



,fe+2 



1=3 



= X\l{jg)(k)w)+0{X''+'), 

where T, e DifP(9A7; "T*M). Given /i = x-'^ {g'g + x^ g'^ H h a;"^;, + 0(x"+i)) as above, according 

to Lemma [ 



2 



2 



!=0 



Hence (3gh e C°°(M;Or*M). Moreover, setting = (-2trgog^)^, 

2l3gh-v° e a;iC°°(M;°T*M)). 
Let = 0;° = -(/( Jg)(0))-it;O - {^trgj^)^, we have 

Jgw(°) + 2pgh - xv^ e x^C°°(M; °T*M). 
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I{jg){k) = 



where = J^tiil^ m 9o]^ + ^^d^trg^g'^)^ . Choose uj^ - -{I{Jg)(l))-'v\ then w^^) - w^^^Hxw' 
satisfies 

Jgw(^) + 2/3g/i e a;2C°°(M;°r*M). 

Since for 1 < fc < n, 

+ nk + 2n 

-fc2 + nfc + 71 + 1 

is invertible, we can repeate this process up to = ti. In other words we can construct a sequence 
of approximate solutions w*^*^-* = uj'^ + xlu^ + ... + x^uj'^ to equation in the sense that 

JgC^W + 2Pgh = x^+\''+^ + 0(a;'=+2), 

where uo^ = uo^^ + ^ and = w^f^ + Note that if fc < n - 1, for < /i < ti 

cjfe = Tfcg^, e C°°(5M) and = Tj^g'^ £ C°^{dM) for some differential operators Tfc with coefficients 
depending only on the components of go and their derivatives up to order k. However for k = n 
then 

Vo=T,,g'o + Togi v^^%,g'^. 
Hence w" = —{I{Jg){n))~^v" satisfies the same property, i.e. 

u;^ = T^g'o + fog'o, Lu^^T^g'o- 
Now = + XLJ^ + ... + x'^uj'' satisfies Jgw(") + 2(3gh = a;"+iw"+i + 0(a;"+2). However, 

r 2n 

/(n+l)=[ 

is degenerate. Notice that {—{xdx)^ + nxdx + n + l){x"-^^ log a;) = — so we can choose 

1+1 _ 1 77"+! — \^ „n+l'^y' 

i=l 

Then cj("+i) = + + (a;"+i loga;)!!; satisfies 

Jg(a)("+i)) + 2(3gh = re {x"+^ log x)C°°(M; °T*M) + a;"+2C°°(M; "T*M). 

Note that for all fc > n+ 1, I{Jg){k) is invertible. Repeating this process again, we can find a refined 
sequence of approximate solutions 

N N 

c^W^^xV+loga; J2 

k=0 k=n+l 

for arbitrarily large with G C°°{dM-,°T* M) and satisfying 

Jgc^t^) + 2pgh = 7^+1 e x^+i^°(M; °T*M). 
By Proposition [4?3l choosing large enough, we have 



Jg"^7^+^ G x"+^C°°(M;°r*Ai') + x ''^'^2'^'" C°°(M;°r*M). 

Then cj = o;'-^-' — Jg"^r^+^ is a true solution to equation (P^ . Moreover, letting 5 = cDoc^a^ + 
Sr=i '^i'iy^j then near the boundary, 

wo = +xul + --- + a;"+ic^o"+i + 0(x"+^)), 

= x-'^ixujl + ■■■ + x"< + a;"+i logxcl^r+i + 0(2;"+'=)), 

where ui^ = ^^trg„g'„. Hence, the dual vector filed X e C^-^iMi^TM) = xC"-'^ (M;TM), which 
implies that X vinishes at the boundary. Hence the 1-parameter family of diffeomorphisms generated 
by X is C"'~^ over AI which restricts to the identity on the boundary. By definition, 

= + C"^-^(M;S^CT*M)) ^ x-^C''-\M;S'^T*M) 

which satisfies the gauge condition Pgh = 0. More explicitly, using the formula 

[S*gU!]al3 = + V/3Wa) = ^ (9„a3/3 + dfjUJa) - T^^W^ 
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and lemmas in the previous subsection, we can compute the components of as 

Tl-l 

fe=i 

ri-l 

fc=i 

n-l 

[5mij = x-\-iol[go\ij + x'Fk + x"{Fn-i - uMij + ^^^oignhj) + 0{x"+')) 

fc=l 

where Fk is a function of a;°, o;^, cj^ for < ^ < n and the components of go and their derivatives. 
Clearly, Fk is linear in for 1 < i < k and < < n. Therefore, writing 

h = x-'^{ho + xhi + ... + + logxff + 0(x"+i)), 

we established: 

(a) ho = g'o- ^(i''so5o)ffo, which imphes that trg^ho = 0; 

(b) For 1 < k < n—1, hk = g'k + Tkg'o for some differential opeator T/. of order k with coefficients 
only depending on go and their derivatives (g^ = if is odd.); 

(c) hn= g'n+ Tng'o + 2o5o- ^0 is of Order with coefficients depending on the components of go 
and 5„; 

(d) H ~ T„+i5q + Ti5q. The coefficients of Ti depends on the components of go,gn and their 
first order derivatives. 

Furthermore, if trg^g'^ = Sg^g'^ = 0, then in the construction of lo, we have — Tk-ig^ for 
< fc < n — 1 and = Tn-^g^ + Tog'^. Hence in this case, (a),(c) and (d) are replaced by 
(a') ho = g'o; 

(c') hr^ = g'„+Tn-ig'o+Tog'o; 
(d') H = Tng'o + T,g'^. 

□ 

5. Linearized Einstein Equation 

After fixing the gauge, the linearized Einstein equation becomes elliptic, i.e. 

Lgh = 0, where Lg = Agh + 2R + 2n. 

In this section, we want to study Lg as was done for Jg and with further analysis of the Poisson 
operator and scattering matrix to characterize the solution to the linearized Einstein equation. For 
this goal, we can study a family of operators instead. Denote 

Lg{X) = Agh + 2R + 2n- X{n- X), for A e 

where Dg is a simply connected neighborhood of [n,oo) in C defined by 

(27) De^{XeC:^X-n>-e, |3?A| < e} 

for some £ > and small. Throughout this section, the parameter A will be in such a domain 
for some s > and small enough. 

5.1. Indicial Operator. For any q e dM, let ([0, e) x U;x,y) be the local coordinates. At the 
boundary, S'^{^T* M)\qJj has a decomposition w.r.t. the metrics g and go, i.e. S'^{^T* M)\qjj = 
y° © © © r^, where 

r° = span{^ (g) 

= {h€ S'^{°T*dM);trg^h = Q), 

= span{n!^(S)^-x-''go}, 

= span{f + : l<z<n}. 
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Note that 'Y^ doesn't depend on the choice of local boundary coordinates {U;y). Hence such 
decomposition is well defined on DM. Moreover, it can be extended to a neiborhood of the boundary 
up to 0{x^), i.e. we can choose A = 2 in Condition [T] That is because 

g x^'^idx'^ + g{x)) = x^'^{dx^ + .go + 0{x^)) 

in a neighborhood of dM . If we replace go by g{x) = go + 0{x'^) in the above definition of Y^, the 
decomposition can be extended to whole neighborhood of dM. The indicial operator corresponding 
to such a decomposition is 



/(L,(A))(,s) 



2n 



2n 

n+1 



Hence for A G D^, the corresponding indicial functions and indicial roots are: 
/o(s,A) = -s2+ns + 2n- A(n-A) =0 

^ So^Aj — 2 ' * ['^) — 2 ' 

/i(s,A) = -s2 + ns + 2n-A(n-A) =0 

^ si(A) =n- A, s\\) = X, 

/2(s,A) = -.s2 + ns + 2n-A(n-A) =0 

^ S2(A)=so(A),52(A)^,sO(A), 

/3(s,A) = -s^ +ns + n+l- X{n- X) = 

_ n-A/(2A-n)^+4n+4 3,.^ _ n+-v/(2 A-ri)^ +4n+4 
^ S3(,Aj — 2 I l-^J — 2 • 

Then if e > is small enough, we can fix a choice of square roots such that Si{X) and s*(A) are 
holomorphic in D,. Note that 5RA = SRs^A) < ^s^iX) < 5Rs°(A) = ^s^iX) < 5RA + 2 for aU A e L>e- 

5.2. Normal Operator. Consider the half-plane model of hyperbolic space (]HI"+^,r), where 



irn+l 



r = x~^{dx^ + dy^). 



At each q G 9M, let ([0, e) x C/; x = y^, y) be a coordinate system. Then S'^{°T*M) hfts to a trivial 



bundle over T+M 



irn+l 



I.e. we can 



use { 



< j < f^} as a basis of 



5'^(°r*lHI"+i). And the decomposition corresponding to 1^° ® is 







^00 


= hii, hoj = /ly = 0, for i,j > 1}, 


#-1 = 


{/. 


E 

i>i 


hi = 0, /loo = hoj = 0, for j > 1}, 


= 


{/. 


^00 


= -nh,i,ho-i = h,j = 0, for i,j > 1}, 




{h 


^00 


= /ly — 0, for «,j > 1}. 


The fibre dimensions di of 


are 






rfo = ^2 = 1, di = — 


n-l) 
2 


+ n 


- 1, ds = n, N = do + di + d2 + d^ = 



(k+1)(k+2) 
2 



Then the normal operator can be expressed as 

n 

Ng{Lg{Xj) = -ixd,f+nxd,-x^Y.dl,+2x 



































E23 





£^31 


£^32 






2n 



2n 

n+1 



- A(n- A), 
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where E,j e DifF^( 



Hofn{W; W^). More explicitly, for h = E"j=o 



a; ' 



da; 



fc=l 



dx 



]=i k=i 

In particular, E{h) = if /i G C°°(IH["+^; #'^) satisfies 6rh — 0. In this case, the normal operator 
becomes diagonal and so effectively reduces to acting on functions. 

Notice that Nq{Lg) = Lr — A.^ + 2Rr + 2n. We deal with the normal operator Nq{Lg{X)) almost 
in the same way as was done for Jg. So we omit the proofs of some statements if they are essentially 
the same as in the last section. 

For any u G a;°ff^(H"+i; C^) for some a € M, let u{x,^) = ^y^^{u){x,^) be the Fourier Trans- 
form w.r.t the variable y. Then u{x, t) is a solution of {Lr — X{n — X))u = if and only if u{x, ^) is 
a solution of the ordinary equation {Lr — X{n — X))u ~ 0, where 

2n 



2n 

n+1 



Lr = -{xdxf + nxdx + a;^|<^p + 



2ixE{i) 



where E is homogeneous of order 1 in ^. Moreover, for any ^ 7^ 0, by changing variable to a; = a;|^| 
and ^ = {Lr — X{n — X))v{x\£,\,£,) = is equivalent to {Lr — X{n — X))v{x,i,) = 0, where 

2n 



2n 

n+1 



-{xdx)'^ + nxdx + x^ + 



2ixE{£). 



Lemma 5.1. For a E s+v^EEl)^ the map 



Nq{Lg)^ Lr:x''H^{W'+'; 
is continuous and positive. 

Proof. For any h G H^{W+\ Ng{Lg){x''h) = a{n - a)x''h + x" Ng{Lg){h). Let h = + h'^ , where 
= ur for some u G C°°(H"+i) and trrh'^ = 0, then Lr{h) = (A^ + 2n)h^ + (A^ - 2)h'^. Hence 

{Ng{Lg){x°'h),x''h)^aL2 = {a{n~ a)h + Lrh,h)L2 

= INh'Wll + llV/i'lli. + {-a^ + na + 2n)\\hXl 

+ (-a^ + na- 2)\\h'^\\]^2 
> 

with equality if and only if h — 0. The continuity is obvious. □ 

Note that since n > 3, [1, n - 1] C "+v^ ) c [ "-v^ , »+v^ ] ^ 

Lemma 5.2. There exist S > and e > 0, such that for a £ {^ ~ ^ + S) and X G D^, where 
is defined by {21^ , the map 

Ng{Lg{X)) ^Lr{X) ■.x"'Hl{W+^-<C^ ®Tl) — ^ a;"L2(H"+i; «) r|) 
is continuous and injective. 
Proof. According to the proof of last Lemma, 



{Ng{Lg{X)){x''h), x''h)^aL2 = II V/i||i2 + (-a^ + na - 2 - A(n - X))\\h\\l2 + 2{n + l)\\h^\ 



,1||2 

r 2 



If (5 > and e > are small enough, then for X e and a G (§ — (5, § + 5), either ^{\ — 5^ - 2 — 
X{n - A)) > or 3(2^ - 5^ _ 2 - A(n - A)) 7^ holds. Therefore the map is injective. □ 
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Let Tg" (M+) be the half density over M+ triviahzed by | jifrl^ • 

Corollary 5.1. Let e and S be the same as in Lemma 1 5. 21 then for a G {^^ — 6,^ + 6) and A G D^, 

the two maps 

Lr{X) : x''Hl{M+; (g) r| (M+)) — > x"L^(M+; r| (R+)) 
are continuous and injective. 

The linear ODE Lr(A)w = has a regular singular point at i = and an irregular singular point 
at a; = oo and no other singular points. So a formal solution at a; = is an actual solution and 
for each formal solution at x = oo, there is an actual solution which has asymptotic expansions at 
X = oo corresponding to it. So to find out all the 2N linear independent solutions for Lr{X) — 0, 
we only need to find out all linear independent formal solutions at a; = 0. Since such solutions can't 
blow up at any point in (0, oo), we can extend them to the half real line (0, oo). 

Notice that for \ = n, the difference between two of the indicial roots may be an integer. Then 
for X G Dg, where £ > is small enough, the 2N linear independent formal solutions, which are true 
solutions, are constructed as follows: 

(1) For the basis e° G = x^ I,^(^\'f{x)e^ and vq — x^ Ki,(^x){^)s^ a-re two independent 
solutions with leading term x'^ ^'^^e and a;*''*-'^-'e respectively, where /^(A)(a;) and Ki,(^x){x) 
are the modified Bessel function of the first kind and Macdonald's function, respectively, of 
order ^(a) = V(2A-»W . 

(2) For any basis G for A G D^, note that /j(s^(A) + fc) = if and only if i = 3, A = n 
and k — 1. So for A ^ n, we can construct a solution = x'^ (-^^e^ + 0{x^ ('^)+^). However, 

blows up at A = n, i.e, is meromorphic in D^, with a pole of first order at A = n. 
Define = d\{{\ — n)u^). Then is a solution with leading term x'^ ^^^e^. 

Similarly for A G D^, /i(si(A) + fc) = if and only ifA = n, i = l,fc = riorA = n, i = 
3, A; = n + 1. For A 7^ n, we want to construct a formal solution ui = a;*^*^^'e^ + 0(a;''i'^''+^). 
Ui is meromorphic in I?^ with a pole of second order at A = n. Define vi = d\{{^ ^ nf'ui). 
Then vi is a solution with leading term x'^'^^^^e^ . 

(3) For the basis G since /i(s2(A) + fc) 7^ for all positive integer fc and all i = 0, 1, 2, 3, 
we can construct two solutions with leading term V2 = x^^^'^'^e^ + 0(a;*2^'*'-*+^) and v"^ = 
^s2(A)g2 _|_ respectively. 

(4) For any basis G first note that /i(s^(A) + fc) 7^ for all fc G N and i 0, 1, 2, 3. Hence 
we can construct a solution = a;* + 0{x). 

For A G -De, /i(s3(A) + fc) = if and only ifA = n, i = l,fc=lorA = n,i=:l,fc=:n+lor 

A = ri, z = 3, fc = n + 2. For A 7^ n, we can construct a solution M3 = x'^^^^^'e^ + 0{x'^^^^^^^). 

U3 is meromorphic in Dg with a pole of third order at A = n. Define V3, — 9|((A — n)^U3). 

Then V3 is a solution with leading term a;*^('^'e'^. 
All the solutions Vi and are holomorphic in for i = 0,1,2,3. On the other hand, at a; = 00, 
the solutions constructed above have leading term either a;^e~^ or x^e~^^. By Corallarv 15.11 the 
behavior of solutions constructed above is 

Vi^x'^e~'^, v^^x^e^ 

as a; — > 00 for i = 0, 1, 2, 3. 

For A G Mn -De, since Lr{X) is self-adjoint w.r.t x~'^~^dx, using the standard method again in 
[CL] we get the Green function for Lr{X), which is 

Gi (A, X, x\ i) = U{\ X, i)V* (A, x', i)H{x' - x) + V{X, x, ^U* (A, x',i)H{x - x'), 

where U (A, x, ^) and V{X, x, f ) are N x N matrices whose columns are linear combination of and 
Vi, namely, 

(28) Lr{X)Gi{X,x,x',i) = x"+^S{x - x'). 
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Moreover, 

(29) C/'(A, X, i)V* (A, x', i) - y (A, x, i)U*{X, x\ = ~\x''-^Id. 

This implies that 

U{X,x,i) = (ao-y",aiw\a2W^, ,03-1;^), V{X,x,£,) = (fco^^o, ^i^i, &2i'2, ^a'l's), 

for some ^ 0, 5^ 7^ for < i < 3. Note U, V, ai, hi can be extended hofomorphicly to 
for e > small and so are U* , V* . Denote U ,V be the extension of C/*, V* . Remark that the 
conjugate relation doesn't holds any more, but only U{X) — {U(X))* y{X) = (y{X))* . Hence Gi can 
be extended to whole such that ([28| and ([29]) still hold in such domain by analytic theorem. 



Lemma 5.3. If e > and S > are small enough, then for all A G and a G — 6, ^ — ■^), 



Gi{X,x,x' induces a bounded integral transformation on x°'L^' 



Lemma 5.4. If e > and 5 > are small enough, then for all X G -De and a G — 5, ^ — |), there 
exists a unique kernel G2{X,x,x' which is bounded on x'^L'^{M.'^]L'^{W^]C^)®Tq{W^)) satisfying 

Lr{X)G2{X, X, x', S) = a;"+M(a; - a;') ® Id. 

Proof. Let G2(A, x, x', ^) = |^|^"Gi(A, a;|^|, ^). The remainder of the proof exactly follows that 
of Lemma 14.51 □ 

Proposition 5.1. If e > and S > are small enough, then for all A G -De and a G — 6,^ — 



Lr{X) : x''H^{W"+^-X^ ® T§ (M+)) — > x" LI{W"+^ ; C'^ < 
is invertible with Green function 

1 



. Tf] (R+ 



)) 



(30) 



G{X,x,x',y,y') 



ier"Gi(A,xicu'ieu>^«(^-^')de 



Moreover, G(A) G ^'o 



-2,.Et(A),-Eb(A) /|nr„_|.l^ , . 



EriX) 



(H"+i; Af(iV, C) ®Tg), where 

sO(A) s\X) + 2 s^{X)+2 s^X) + 2 

.sO(A) + 2 si (A) s2(A) + l s^iX) + l 

sO(A) + 2 si(A) + l s^{X) s\X) + l 

sO(A) + 2 si(A) + l s2(A) + l .s^{X) 



and Eb{X) = Et{XY 



As in the last section, let x be the boundary defining function of the hyperbolic ball -6"+^ and p 
and p' be boundary defining functions for Fq. Let ei(A) — (s'(A) + 1, s'(A), s*(A) + 1) be the 
constant vector corresponding to the i-th column of Et{X). 

Corollary 5.2. For X G D^, where e > is small enough, the map 

3 



G(A) : G°°(B"+';C^(^r|) 



is continuous. 

Corollary 5.3. For X G D^, where e > is small enough, the map 
G(A) :G°°(TV;M(7V,C)0r|) 

is continuous. 
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5.3. Parametrix and Resolvent. 

Proposition 5.2. If e > and S > are small enough, then for all \ <E and a G — (5, ^ — |), 

Lg{\) =Lg- X{n - A) : x^HUM- S\^T*M) ® tJ) — ^ a;''L^(Af; S^CT*M) (g) T^) 
is Fredholm. There exist right and left parametrices 

such that 

Lg{\)Qr{X) -Id = Rre ^--'°0'^«(^)(M; S\°T*M) ® rj), 

Qii\)LgiX)-Id = Ri e %°°'^^'^^^'°^iM;S^°T*M)®rl). 
QiiX),Qr{X), Ri{X), Rr{X) are holomorphic in Di^. 

Proof. Based on the analysis for the normal operator in the last subsection the existence of the right 
parametrix Qr{X) follows directly by Theorem 13. II By Proposition 13.11 

Qr{X) : x^LliM; S^{°T*M) ® r|) — » x^H^iM; S^{°T*M) (g) r|) 

is bounded and 

Rr{X) = Lg{X)Qr{X) -Ide ^--'-.^B(^).(Af ; S\°T*M) ® rj). 
Moreover, since Lg{X)* = Lg{X), we can define Qi{X) = Qr{X)* . Then 

RiiX) = Qi{X)LgiX) -Id= (Lg(A)Q,(A))* -Ide m,---^^-(^).-(M; ^^("T^M) ® r|). 
Note that Rr{X) and i?;(A) are two compact operators on x°'Lq{M; S'^{'^T* M)(^Tq). That is because 

Ran(i?^(A)|^„i2) c {IT; S^^T* M) (g)Tl), 

1=3 

Ran(i?i(A)L.i.) C ^ ^^("r^M) ® Tj), 

i=Q 

where ei(A) is the i-th column of Et{X), and the inclusion maps 

s/^{M;S'^('T*M)®tI) ^ x''Ll{M;S^{°T*M)®Tl), 

s>/^^(^) (M; S^i^T* M) ®tI) ^ x" LI{M; S^{°T* M) ®tI) 

are compact for X <E and ae(^ — (5, ^ — |), where e > and (5 > are small enough. Hence 

Id + Rj. and Id + Ri are Fredhom operators on x'^LliM ; S'^{^T*M) ®T^)^ which implies that Lg{X) 
is Fredholm. □ 

Lemma 5.5. If e > and S > are small enough, then there exists Aq G M n ZJ^ such that for all 
A e M, A > Ao or A e De, 3A 7^ 0, and for a S (§ - 5, f - f ), 

Lg{X) : x"Hl{M- S^CT*M) (g) fJ) — > x^LKM; S^{°T*M) ® fJ) 

is invertible. Moreover, Lg{X)-^ G %^''^^^^^'^''^^\M; S'^{^T* M) ®Tl). 

Proof. Since dM is compact, we can choose a finite coordinate covering {Ui : i = for 
dM. Hence we get a finite number of coordinate patches {[0, e) x Ui : i = l,...,p} which cover 
[0, e) X dM for some e > and small. Since M — [0, e/2) x dM is compact, it also can be covered 
by finite coordinate patches. Thus we have fixed a finite coordinate covering for M. By directly 
computation, there exist some constant C such that any sectional curvature < Cx'^ in any 

coordinate patch. So for any h € H^iM; S^{°T*M) ® f|), 

\\Rh\\<C,\\h\\L. 
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for some constant Ci. Choose Ao = "+^"^+^'^^ . Then for any h e H§{M; S^{°T*M) (g> T§) and 
A e [Ao,oo) 

{LgiXo)Kh)Li = \\Vh\\l, + {2Rh,h)Li + i2n-X{n-X))\\h\\L2 
> l|V/i||i2 + (2n + A(A - n) - 2Ci)||/i||i2 > 

and equaUty holds if and only ii h — 0. Hence Lg{X) is positive and thus invertible. If A € I?e and 
SJA 7^ 0, then 9(2n — X{n — A)) 7^ 0. Hence Lg{X) is injective and surjective since its adjoint Lg{X) 
is injective. The second statement is shown by a proof similar to that of Proposition 14.41 □ 

Proposition 5.3. If e > and S > are small enough, then for all X E and a E — S, ^ + 6), 
the resolvent Rg{X) = (ig(A))-i E %^'^^^^^''^''^^\M; S^{°T* M) (g)rl) and as a map, 

Rg{X) ■.x''Hl{M;S'^{^T*M)®Tl) — > a;°Lg(A/; 52("T*M) (g) r|) 

is meromorphic in with poles of finite rank at each point of {A E : X[n — A) G o'pp(Lg)}. 
Moreover, {A £ -D^ : A(n— A) E app{Lg)} is discrete and finite. If Xi E Dg, and Xi{n — Xi) E app{Lg), 
then the residue of Rg{X) at Ai is 



mx,) ^ ±-Jjgix)d. 



with kernel k(H(Ai)) = X]^=i Wi(-^i) ^"^(Ai) where Uj{Xi) E X]^=o ■'^'^'^^^^ '^^^ Lf^-engenvectors 
of Lg and 70 is a small circle around Ai. 

Proof. Let Aq be the same as in the last Lemma. A similar proof as in Proposition 14.41 shows that 
for A E D^, 3fiA > Ao or A G D^, 9A 7^ 0, 

RgiX) E %^'^^'^^'>'^"^^\M;S^'^T*M)(^tI) 

is holomorphic. By Analytic Fredholm Theorem, Rg{X) — (Lg(A))~^ is meoromorphic on D^. Hence 
there are only finite poles in D^. By choosing a limit from upper half plane, we have 

Rg{X) E %^'^^''^^'^^^^\M;S^''T * M)(g,Tl) 

for all A G -Dj, Ai(n — Ai) ^ app{Lg). The finite rank property of poles follows directly from the 
Fredholm property of Lg(X). □ 

Corollary 5.4. If X E with e > 0, small and G N satisfies N > 2n, denoting by ei{X) the i-th 
column of Et{X), then 

3 

i?g(A) : x^C^{M;S'^{^T*M)®tI{M)) — > ^ ^'^•(^^(M; 52(°T*M) ® fJ (M)) 

2=0 

is bounded and meromorphic with poles of finite rank at {A G : X{n — A) G (Tpp(Lg)}. 

Proof. For N > 2n, a;^C°°(M; S^{°T*M) (g) r|(M)) C x^H^iM; S^i°T*M) (g) F|(M)) where a E 
(f — 5, f + (5) with S small. Hence if {X ^ : X{n — A) G crpp{Lg)}, the map 

Rg{X) : x^C°°(M;S^CT*M)(E)TI{M)) — > x^LKM; S^^T* M) (g)Tl (M)) 

is well defined and continuous. Moreover, for any / G C°°(M; S'^i°T*M) (g) f|(M)), view Rg{X) as 
its Schwartz kernel, then 

RgiXjf^irURaWi^y)) 
where ttl and nji are the projection map from M Xq M to the first and second copy of M. It follows 
immediately that the image is contained in XiLo ■b^'''^^\M; S^{^T*M) (g) (M)). □ 
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5.4. Poisson Operator. Let > 2n be a fixed integer in this subsection. Recall = {X E C : 
3?A > n — e, |SA| < e}. Fix e > and small enough such that all the lemmas and propositions in 
Subsection 5.1-5.3 hold. Then for A € D^, define maps 

$.(A) : C°°{dM; f ® T^ipM)) — > S^{°T*M) (g) r| (M)) 

for i = 0, 1, 2, 3 by formal series construction, as we do for finding formal solutions for the ODE in 
Subsection 5.2, such that 

Lg(A)$i(A) : C°°{dM;y' (g)r^dM)) — > x^C°°(M; S'^(°r*Af) «) r| (M)) 

is continuous. Namely, for any ho £ C°°(5M;1^* <g T^dM)), 

N 3 
k=l 1=0 

where uu € C°° {DM; S'^{°T*M)\gjj (g) T^dM)) for aU kj. In fact we can construct this formal 
series up to arbitrary large order. However, we only need finitely many terms here due to Corallary 
15.41 Now, for any A € D^, we can define the Poisson operator following Graham and Zworski in 
[GZ], i.e. for i = 0,1,2,3 

F,(A) = $,(A) - Rg{X)Lg{X)<i>,iX). 
Hence Lg{X)Pi{X) ~ 0. The mapping property of Pi (A) follows immediately. 

Lemma 5.6. For i = 0, 1, 2, 3 and X G D^, denoting by ei{X) the i-th column of Et{X), then 
P,(A) ■.C°^'{dM);'r' (gr^dM)) — ^ a;"'(^V°(M ; S'2(°r*M) (g) (M)) 

3 

+ ^/'^•(^)(M; S^CT*M) (g) r| (M)) 

i=0 

is continuous and meromorphic with poles of finite rank at {A € I?£ : A(n — A) € app{Lg)}. Namely, 
for any ho G C°°{dM); f (g) T^dM))), 

3 

3=0 

where F, Go, d, G2, G3 £ j^°(M; S^{°T* M)) (g vl (M)) satisfying F\J^\-^gjj = ho. 

Notice that for A G D^, iam{^s°{X),^s'^{X),^s^{X),^s^iX)} = 5RsHA) = 5RA. Moreover, if 
A ^ 71 + ^ , we have 

{|,s,(A) - ,s,(A)|, |s'(A) - s,(A)| : z, J = 0, 1, 2, 3} n N = 0, 

which implies that no logarithmic terms have involved in the construction of formal solutions to the 
ODE and hence we can choose X = in the definition of p?]) . i.e, we can replace all the -^phgiM) 
by G°°(M) in the previous part of this section. Namely, 

Lemma 5.7. For i = 0, 1, 2, 3, A e and X ^ n + the two maps 

3 

Rg{X) : x^C^iM; S\°T*M) (g rj(M)) — > ^ x"^ C°° (M ; S\°T*M) g) r|(Af)) 

i=0 

Pi{X) : C°°{dM; (gV^dM)) — > cc"' G°° (M; S^^T* M) (gvl (M)) 

3 

+ ^a;"'(^)G°°(M; S\°T*M) (gvl (M)) 
j=o 

are continuous and meromorphic with poles of finite rank at {A G Df, : X <^ n + A(n — A) G 
app{Lg)}. Hence for any ho G C°° {dM]Y^ g) T^{dMj), there exists a constant < c(A) < 1 such 
that 

Pi(A)/io = x'^-^F + x^G + 0{x^+^^^^) 



dx ^ 

.n+l 
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where F,G e C°°(M; S^{°T*M) O r|(M)) with 

Proposition 5.4. For X G D^, X ^ n + ^ and X{n — X) ^ app{Lg), the map 

Pi(A) : C°°{dM;y^ (E)r^dM)) — ^ a;"-^C°°(M; 52(°T*M) (g) rj (M)) 

3 ^ 

+ x'''^^^C°° (M; S'^CT*M) O r| (M)) 

3=0 

has kernel k{Pi{X)) = Rg{X)x'~^~^^\dx'\~^x'=o- a distribution on M Xq dM with 

singsupp(«;(Pi(A))) C dM Xq dM. 

Proof. For any ho G C°°(OM; O Ti{dM)), f e x^C°°(M; S'^{^T*M O r|(M)), by last lemma, 

ui = Pi{X){ho) = + x^Gi + 0(x^+^(^)), W2 = Rg{X){f) = x^G2 + 0(a:^+^(^)) 

for some < c(A) < 1 and Fi,Gi,G2 € C°°(M; S'2(0T*M) (g)r|(M)). Then for X G RnD^, 
A ^ n + ^ and A(n - A) ^ app{Lg), 

/ Pi(A)/lo/ = / UiAgU2 - / AgMiU2 

= lim / Mi(Vu2-'^)— / (Vui • n)u2 

= (2A-n)/ PiG^^c^+^xr^laM 
7aM 



= (2A-n) / /ioa;-^+'^i?g(A)/|da;|-5|aM 

JdM 

= (2A-n) / (i?g(A)a;-^+'^|da;|-5)/io(g)<5o(a;)/ 
Jm 



The result extend to whole domain {A € : A ^ n + A(n — A) ^ app{Lg)} by analytic 
continuation. Since -Rg(A) only has singularities on d,iagf){M) and cornormal singularities on the 
boundary surfaces of M Xg M, k(Pi(A)) as a restriction on one boundary surfaces must be smooth 
when a; > 0. □ 

If Ao = n + I G for some / G No, then for any /iq G C^{dM; O T^dM)), 

Pi{ho) = x'^-^" F + x^" {\ogx)Hi + x^°+\\ogx)^ H2 + 0{x^°) 

where F, iJi, iJa € C°°(M; S^{°T*M) ® (M)), /fa = if Z 7^ and 

J^l^r^laM = /io, H,\^\-i\^jj G G^{dM;r' ® T^dM)), 

H2\^\-%-MGG'^{dM-r^®T-2{dM)). 

^\dM H2\-p^\ ^\g-M determined by the formal construction of $1 up to order 
N, i.e. there exist some differential operators 

An+i GD-^+\dM;End{y^)(^T'^{dM)), An+iGT)i^+^+'^{dM-Hom{'r^,y^)®T'^{dM)) 

such that Hi\J^\-^\gj^ = An+iho, H2\J^\-^\gM = An+i+ih^. 

If Ao G Dg, Ao(n — Ao) G app{Lg), then the Poisson operator has residue of a finite rank operator 

at Ao 

Resx=Xo{PiW) = (2Ao - n)n(Ao)a;'-^U/=o. 
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5.5. Scattering Matrix. For A G D^, \ ^ n + ^ and \{n — A) ^ (7pp{Lg), we can define the 
Scattering matrix 

Si{X) : C°°{dM] ® Y^{dM)) — > C°^{dM- 1/^ ® Y^{dM)) 

by 

S'i(A)/io = Gx'^\dx\~^Qj;j. 

Since the choice of smooth boundary defining function for dM remains free, it is more natural to 
view 5*1 (A) as a map 

S'i(A) : C°°{dM\ ® Y^{dM) ® \N*dM\^-^) — > C°^{dM- ® Y^{dM) ® \N*dM\^) 
Proposition 5.5. For e > 0, small, the scattering matrix Si{X) is meromorphic in with poles 



ai n + ^ U {A e : A(n — A) ^ app{Lg)}. At any pole Aq, we have 



' 7r(Ao) + forso = n+5, 



where «;(^(Ao)) = (2Ao - n)a;-^"n(Ao)x'-^° | J^^^t^^ |^ | 



dMxdM- 



To find the kernel of 5'i(A), for A G D^, X ^ n + and A(n — A) ^ app{Lg), we study the 
asymptotic behavior of k{Pi) first. Recall the blow down map 



P : M XqOM — > M X dM, (3 : dM Xq dM — > dM x dM. 

The front face for /3 is n and the front face of /3 is F n i? n T. For any {q, q) E diag{dM) and a 
neighborhood U C M x dM around it, we can choose local coordinates in C M Xq dM as 

R^\x' + \y-y'\hl Y = y^y', ^ = pp = {1 + HY'^ ■ 

By Proposition O for any ij = 0,1,2,3, /?*(k(Pi(A))) ^ R-^Fv^vq, where 



F 



Coo(A)p^°(^)Foo + Coi(A)p^'(^)+2Foi + Co2(A)p^'(^)+2Fo2 + Co3(A)p^'(^)+2Fo3 

Cio(A)p^''(^)+2Fio + Cii(A)p«'(^)Fn + Ci2(A)p^'(^)+iFi2 + Ci^{X)p'"^^^+^Fi^ 

C2o(A)p^''(^)+2Fio + C2l(A)p^'(^) + lF2i + C22(A)p^'(^)F22 + C23(A)p^'(^) + lF23 

C3o(A)p^''(^)+2i^30 + C3l(A)p^'(^) + lF3i + C32(A)p^'(^) + lF32 + C33(A)p^'(^)F33 



with Fij = Fij{p, R,uj) e C°^{AI Xq dM) and the coefficients Cij{X) are meromorphic functions of 
A. Let Sij{X) — s-'(A) + 1 — Sij and 

T,j (A) = x^'^^ (^)/3, {R-^F,j)v ® I/O = /?* {R-^-"^' (^^Fy ® i/q. 

Since Sy > A except for i = j = 1, to study the scattering matrix Si (A), we only need to consider 
the asymptotic behavior of Ty(A). 

Proposition 5.6. For X e Z3e, A ^ n + A(n — A) ^ app{Lg) and for i,j = 0, 1, 2, 3, Ty(A) is a 
distributional section ofT'i{M x dM), which has conormal singularity at diag{dM). Moreover, it 
has asymptotic expansion in x as x I in the sense that for any f € C°°{dM x dM), 



dx 



where Hx,Gx G C°°{[0,e)) depend holomorphicly on X. Moreover if i = j ~1, 

Hx{0) ^ {(3*{R^'^^F^j\p^o)vo (g) Vo,fvo ® l^o)gMxOM 

G\iO) = {{x'^^^^nj\pnB,PF^'~W)FnBSoiY)i^o(^i'oJ''o(^i'o)oMxd-M 
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Proof. The proof follows step by step from the proof of Proposition 4.1 in [JS] by Joshi and Sa 
Barreto. We only need to do for such / that supp{f) is contained in a small neighborhood of 
(q, g) G diagipM). In the radial coordinates around {q,q), 

RdR = P*{xdx + Y -dy), d^R^p, xd^p^ p{l- p^). 

For any m S No, 

KUi^d, ~ k){xd, + YdY + \ + s,,{X) - k)n,{X) = 
where F^^ = R, uj) eC°°(M Xq dM). Let 

<^)=l_l MR'^'"''^^^F,j){x,v,Y)f{y,Y)dYdy, x>G. 

JdM JR" 

Using the identity div(Yg{Y)) = ng{Y) + Y ■ dygiY), we deduce that 

nfcLola;^:, -k){xd^-n + \ + .s,y(A) - k)u{x) < C\xr+\ x > 0. 

Since A + Sij{X) is not an integer, the ODE theory shows that for any large m, there exists some 
constant m' independent of m and integers Pi , P2 , such that 



Pl P2 

k=0 1=0 



u{x) -J2^''dk - 



< C\x 



where dk,d[ are constants. Hence H\{Q) — do, Gx{0) — d'g. 

Furthermore, if i = = 1, then (A) — A, which can be extended holomorphically over C. Hence 
Til = /3* (R~'^^Fii)i^ (g) lyQ can be extended meromorphically over C. As in [JS], the formula of G\{0) 
and H\{0) follows directly by analytic continuation. □ 

Corollary 5.5. For A € D,,, the kernel of the scattering matrix S'i(A) is a meromorphic family of 
pseudodifferential operators of order 2A — n, whose kernel can be expressed as 

KiSi{X))^x-^+'^Rg{X)x'-^+'-^\dxdx'r-^U^,,^o^\Y\-^^F{y,Y) 

for some smooth function F £ C°°(9MxR"), and with poles at n+^n{A G : A(n— A) G app{Lg)} 
of residue as in Provosition \5.5l The principal symbol of Si{X) is 

a2A-„(5i(A)) - ^y^^(|yr2>F(A,y,r/|r|,o)), 

where F{X,y,Y/\Yl\Y\) = p-^Rg{X)p'-X=p,^o and F{X,y,Y/\Y\,0) = p-^Rg{X)p'-^\R=p^p,^o- 

Proof. For A G D^, X ^ n + ^ and X{n — A) ^ a-pp{Lg), 5Rsy (A) > 5ftA + c for some constant c > 
except for i — j = I. Hence the scattering matrix >S'i(A) only involves rii(A). The formula for 
the kernel follows directly from the preceeding Proposition. Hence «;(iS'i(A)) only has a conormal 
singularity at diagidM). Near diag{dM) 

^iS^{X)) = \Y\-^'F{X,y,Y/\Y\,\Y\). 

Therefore Si (A) is a meromorphic family of pseudodifferntial operators on dM of order 2X — n for 
XeDe. □ 

Corollary 5.6. //S = {ho G C°°(9M; r|) : Sg^ho = 0}, then for X G D,, S'i(A)|s is elliptic 
and meromorphic in with poles n + ^ U {A G -Dj : A(n — A) G crpp{Lg)}. In particular, in this 
case, Xq — n is not a pole if Lg is invertible. 

Proof. Let E — Pi(S), then the normal operator of Lg{X)\^ is diagonal. Hence the Green function 
for Nq{Lg{X)\^) is a function of {R,p,p') by standard result for Ag acting on functions as in [MM]. 
By [JS] and the preceding Corallary, for A G D^, X ^ n + No and A(n — A) ^ app{Lg)}, 

r(A-i) 

□ 



0'2A-n('5'l(A)|s) = 2" ^^-p^ —a2\-n{'^go 
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h = 



6. Proof of Theorem 11.11 

Proof. Under the hypothesis of Theorem ll.il given any family of boundary metrics go{s) G Met°° , 
for s S {—9,6), 9 > and small, there exist a family of C^'" compact conformally Poincare- 
Einstein metric g+{s) such that gois) € [x'^9+is)\rpgjj] and g+{0) = g is smooth. By Proposition 
2.1i there exist a family of C^'" diffeomorphisms $s over M, which are C^'" in Af , such that 
near the boundary 

~g+is) = <^>:g+{s) = x-'{dx' + gois) + x^g2{s)) + ■■■ + x^g,,{s) + 0(x"+i). 
Moreover, $0 is smooth. Hence g = .g+(0) is smooth. Then 

dg+{s) 

dx 

By Proposition 14.51 there exist a family of C"~^ diffeomorphisms ^'^ such that = Id, '^s\qm = 
/c?g;jg and h = ^|s=o(^s5+(s)) satisfies /3gft- — and Lg/i = 0, where 

h = a::"^(/io + x^h2 -\ h a;"/i„ + 0(x"+^ logx)); 

here ho — g'g — ^i?'go5o ^'^'^ ~ 9'n ~ '^iffo '^'^'^ some differential operators Tn of order n with 
coefficients depending on the components of go and gn and their derivatives. Moreover, if trg^g'^ = 
Sg^g'o — 0, then T„ reduces to a differential operator T„_i of order n — 1. 

Recall that = {A S C : 5ftA > n — e, |3A| < e}. The uniqueness assumption implies that Lg 
is invertible. Hence for e > and small, the resolvent family of Rg{X) and Pi (A) are analytic in 
a neighborhood of A = n in D^. According to the decomposition of S^{'^T*M)\qjj in Section [51 
x-^ho G C°°(aM;ri). Hence 

hi/ — Pl{Tl){x^'^holy'o)- 

Recall that v = \dvolg\^ , = \dvolgo\^ . Note that for A e C, < |A — ri| < e, where e > is small, 

Pi{\){x-^hoVo) - + x^G + 0(a;^+") 

for some c > 0, where F\^\-%-j^ = x'^hovo, G\^\-%-j^ G C°°(9M; ® P^ (SM)). The 
fixed half density sections also satisfy i'\ Ialv7 ~ ^o- Hence we can omit them and only write 

h = Pi{X)ix-^ho) = x''-^F + x^G + 0(x^+'=) 

where F,G e C°° {dJl; S^{°T* M)) and F\gjj = x^^ho, G\gjj = Si{X){x'^ho). By Proposition 
15.51 i?esA=n>S'i(A) can only be some differential operator which corresponds to the x" logx term in 
h. However, by Proposition 14.51 after fixing the gauge, the largest logarithmic term is not bigger 
than a;"+^logx, i.e. Si{X){ho) must have no pole at A = n for all ho G x'^C°° (dMiT^). Hence 
Si{n) : ho hn is a pseudodifferential operator of order n. Therefore, 

g'„ = hn + Tgg'o = Si{n){go - ^trg^^g'^) + Tng'o = ^n5o 

where T„ is a pseudodifferntial operator of order n, which is exactly the linearization of the Dirichlet- 
to-Neumann map DN' . 

Furthermore, let S = {7 G C^{dM;T*dM) : trg.j = Sg,j ^ 0}. Ug'o G S, g'^ = Si{n)g'^+T^^ig'a 
for some differential operator T„_i of order n— 1. By Corollarv l5.61 S'i(n)|s is elliptic. Hence DN'\^ 
is elliptic with the same principal symbol as S'i(n)|E. 

If g'o is a general symmetric 2-tensor on dM, we first modify g'^ to be an element of E by a 
linearization of a family of diffeomorphisms and conformal changes, i.e., we want to find ho G S, 
u G C°^[dM) and w G C°°{dM;T*dM) such that 

f ho ^ g'o + ugo + S*g^uj 

(31) <^ trgjio - 

[ Sg^ho = 

From the first two equations, we get u = ^(^^goSo "^''so'^go^)- Replacing u by ^{trgggQ — trggSg^^uj) in 
the first equation and using the third one reduces the problem to finding a solution to the equation 

(32) Sg^tiS^u; = Sg„tfg',; 
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here 'tf means the trace-free part of a symmetric 2-tensor. Let A = Sg„tfS*^Lu. Then 

A : H^{dM;T*dM) — > L^{dM;T* dM). 

is continuous. It's easy to check that A is an elhptic operator of second order. For any q £ dM, 
choose normal coordinate neighborhood {U;y^, y") around q. Then the principal symbol of A at 
q is 



For = (1,0, ...,0), a{A){^o)\q is invertible. Hence a{A){£_)\q is invertible for all ^ ^ since it's 
invariant under the action by 0{n). In fact 



a2iA)-'io = m- 



V 



n-2 
n-1 



6Ci 



Furthermore, A is self-adjoint. Hence Ran{A) = Null{A)-^. Moreover, for any uj e H^{dM,T*dM), 

{ALj,cj)L2^\\tfSlu;\\i,>0. 
This implies that for any oj E Null{A), ~ 0, so 

Therefore, SggtigQ G Ran{A). Furthermore, by standard elliptic theory for compact manifolds, we 
can construe a generalized inverse B, which is an elliptic pseudodifferential operator of order —2 
with principal symbol a-2{B) — (y2{A)^^ such that 

AB ^ Id- H, 

where H is the projection map onto Null{A), which is obviously of finite rank due to the Fredholm 
property. Let u — BSggtfgQ, then it is a solution to (|32|) . Of course, it may not be unique. However, 
we only need one. Using this solution 

ho^ 9o- ^{trgog'o - trgJl^Bdg^tfg'o) - S*g^B6g„tfg'o = ifg'^ - tfSl^BSg^tfg'o := Qg^g'^ 

where Qg^ defined above is a pseudodifferential operator of order 0. The linearizaion of properties 
^ and ([6]) for the Dirichlet-to-Neumann map shows that 

dAfig'o) = dAf{ho) + ^tig-a + S*uj 



□ 



Therefore, an{dU) = ^kY^cmi^lo^go)- Theor em il.ll is proved. 

Finally, we check that the principal symbol formula above is exactly a generalization of Graham's 
result. For n > 5, odd, let 

W : Met°"{dM) — > C°" {dM ; (g)^T* dM) 

denote the Weyl tensor and 

W : C°°{dM; S^{T*dM)) — * C°°{dM;(S)'^{T*dM)) 
its linearization at go. By direct computation 

a,{W*W) = ^a,{Aleg„). 

For n — 3, let 

C : Met°°{dM) — > C°"(dM; S^{T*dM)) 
be the Cotton- York tensor and 

C : C°^{dM- S^{T*dM)) — > C°^{dM- S^[T*dM)) 
be its linearization at g^. Then C is self adjoint and 

cT3(|c|) = a3(Aiegj, 
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where |C| = {C*C)^ . Recall the result by Graham in [Grl] 

Theorem 6.1. Let M = 5"+^ and dM = 5". Then the linearization of the Dirichlet-to- Neumann 
map at the hyperbolic metric go — d9^ is 



(33) dj\f = \ «^*^('^5o + ci) • • • (Ago + c,„)^Ago +c„+i, ifn>5, odd; 

[ 3 |C|, if n — 3, 

where m = ^^y^, a 7^ and q > /or 1 < i < rn + 1 . 

It is now clear that our result is compatible with Graham's. Moreover, we can compute the 
constant a in his formula, namely, 

_ (n~2)r(-f) 

" 2«(n-3)r(f)' 
7. N-PoiNT Problem 

The Dirichlet-to-Neumann map is of high nonlinearity provided that it is well defined. In the 
previous sections, we investigate its linearization as one way to describe properties of the map. In 
this section, we continue to the A^-point problem, to give a better description. 

To set up the problem, suppose the Dirichlet-to-Neumann map is well defined in a neighborhood 
of go G Met°°{dM). Let O be a small neighborhood of in the parameter space and g+{s) E 
C°° {O; PE^'") a smooth family of C^'" conformally compact Poincare-Einstein metrice. g = g+{0) 
is smooth. As in Section [21 via a smooth family of diffcomorphisms, we can assume 

g+{s) - x-\dx^ + go{s) + x^g2{s) + ■■■ + x"gn{s) + 0(x"+i)). 

where g^{s) € C°°{0;C°°{dM;S^{T*dM))). For simplicity as before, let 

g = .g+(0) - x-\dx^ + g„ + . . . + a;"g„ + 0{x"+^)) 

be a fixed smooth metric. Let j3 = (/3i, (3n) with (3i e No and denote 

g^^Hs) = d^g+is), g^'\s) = d^gois), gifHs) = 5f g„(s). 

By assumption, gn{s) ~ DN{go{s)) is well defined. Then clearly, g'lf\o) is determined by gj^\o) for 
all < 7 < /3. From previous sections, gif^O) linearly depends on glf\o) via a pseudodifferential 
operator same as the first order linearization. In this section, the dependence on g^^"^ (0) for < 7 < 
is studied. For simplicity and as a start of induction, we only do the two point problem, i.e. |/3| — 2. 

In this case, we can assume N = 2 and WLOG, let (3 — (1, 1) in the two point problem. Consider 
the second order linearization of the Einstein equation, 

(34) lJ'^^\Q) ~ 2j;/3g.g(i'i)(0) + i?(.g(i-")(0), g(°^i)(0)) = 

where i?((7^^''^\ g^'^'-'^') is a symmetric quadratic term, which is a differential operator of order 2 with 
respect to each variable. For any p £ M , choose normal coordinates (W; y^) around p, then for any 
h,k£ C°°{M;S'^{°T*M)), 

E^,{h, k) = -h'^'id^T^^eik) - d,T^,g{k)) - k'^'id^T^Mh) - d^T^^h)) 

-g'^P g'l'^ [dcrhp^T ^^^0{k) + dakp^T ^^^e{h) - d^hp^Vp^ae{k) - d^kp^V ^aoQi)) 

+ |9"V'(rM.p(/i)r.e7(fc) + T^,p{k)T„e-,{h) 

-Tf^a-fih)T^ep{k) - Tpcr'f{k)Ti,gp{h)) 

at p. Rence E{g^^'°\g^°^^'>) e {M; S^{°T* M)). 

Say a quadratic form Q{h, k) satisfies the asymptotic property, if for some e > 0, writing 

h= {ho + xhi + --- + a;"/i„ + 0(0;"+^)), 

k^{ko+xki + ---+ x"kn + 0(:e"+^)), 
where h, h, e C°°{dM; S^{°T*M)) for < i < n, then 

g(/i, k)^{Qo + xQi + --- + x"g„ + o(x"+^)) 
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with Qi e C°°{dM; S'^i^T* M)) for each < i < 71, which is a sum of quadratic terms in {hj, ki) as 
a differenial operator or pseudodifferential operator of order k with respect to each variable. Denote 
ord((5) = k. Obiously, E satisfies the asymptotic property with ovd{E) — 2. 

Equation is not eUiptic with respect to ^'^^'^^(O). To make it eUiptic, we have to make a 
suitable choice of the gauge. Following the ?? in Section [4l there exist a smooth family of vector 
field y(si), with dual 1-form ro(si) for (si,0) € O, which generate one parameter families of 
diffeomorphisms for each si, such that if letting 5(51,52) = '^l^g+isi, 82)^ then g(°'^)(si,0) = 
5('''i)(si,0) + ^*(^^ 0)n7(si) satisfies 

(35) /?g(.i,o)5^°''nsi,0) = 0. 

Note that g(si,0) = 5(si,0) for aU (si,0) e O since %^ = Id. Hence g^^^°\Q) = g'^^'°\Q). Similarly, 
for each S2, (0, S2) G O, we can choose a smooth family of vector fields X{s2), with dual 1-form 
u;(s2)j which generate one parameter familes of diffeomorphisms <I>jJ for each S2 small, such that if 
letting G(si,S2) = m~g(,si,S2), then G(1'°) (0, S2) = 5'''°H0, S2) + ix(.2)5(0, S2) satisfies 

(36) %o,..)G(i'°)(0,s2)=0. 

Note that <i>Q^ — Id, so G(0, S2) = g(0, S2) for all (0, S2) G O. In fact, it is only necessary that 137(51) 
and w(s2) satisfy 

(37) (A<,(,,,o + ")t^(si) + /3s(si,o)3^°''Hsi,0) = 0, 

(38) (Ag(o,,,)+n)cc>(s2)+/33(o,..)ff<''"H0,S2) - 0. 
By direct computation at s = 0, 

G(°^i) = 5(°'i)-<5;(A, + n)-i/3,g(0'i), 
Gd.o) = 5(1-0) -<5;(A, + n)-V,5(^^°\ 
G(i-i) = 5(^'^^-'5;(A, + n)-V.5(^'^)+g(.9("\5(^'°)), 

/33(G("-i)) - /33(G(i'0)) = 0, 

/^^(Gfi'i)) = i?(G(°-i),G(i'")), 

where Q and i? are quadratic terms satisfying the asymptotic property with ovd{Q) = 2, ord(i?) = 1. 
Moreover, 

g(^(04)^ g(i,o)) = (Qo + xQi + • • • + x'^Qn + 0(a;"+^)) 

with g, = a;-2^^.^fc<,Gjfe(y)V\g^°'^^V^g^^'°) G G°°(9M, S'2(0T*M)) for each Q <i <n~l. In 
particular. 



= ^ G,,V'=5^°'^)V^5^^'°) + {To{go,gn,gngir' + (To(5o,ff^°^))ffi°'^) 

j+fc<n-l 



where Tq and To are differential operators of order with coefficients depending on the components 
of 50) 5rii So"'^'' (resp. ffo: 5n, ^0^'°'') and their derivatives, which are linear in ^q"'^'' (resp. g^^'^^)- 

Since G(si,S2) is again a smooth family of Poincare-Einstein metrics, G^^'^^ G*^i'°\ G^^'^' also 
satisfy the equation Hence 

Lg(G(i^i))-^(G(i-"\G("^i)) = 

where £;(G(i^°),G(''^i)) = -E{G^^'°\G^°'^'>) + 26*gRiG^^'°\G^^'^'>) also satisfies the asymptotic prop- 
erty with ord(i?) = 2. Moreover, with respect to the decomposition of S^{^T* M)\gjj in SectionEJ 
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UiEo = 0. Writing H = G^^'i), h = G^^'") and k = G(°'1\ 

= Ha + xHi + ■■■ + x'^H„ + 0{x''+'') with 

/i = ho+xhi^ h x"/i„ + 0(a:"+') with 

ho=x-hf,,gi'^'\ K^x-'{9(^''^~T„gl'% 
k = ko + xki-\ h a;"fc„ + 0(a;"+^) with 

fco = x^'^Ug„g^°'^\ fc„ = a;"2(^(04) _ T„5o'^), 

where T„ is a differential operator of order n with coefficients depending on the components of go 
and gn and their derivatives. By Section [H h = Pi{ho) and A: — Pi{ko), where Pi is the Poisson 
operator defined in Section [5l Hence Pi is defined by its kernel 

p{m;y') G C-°^'(M x dJl; Hom{y\ S^{°T* M))). 

Note that we omit involving the 0-half density here since the backgound metric g and g^ trivialize 
Fq (M) and (9Af ) respectivley. Denote m = (a;, y) in a neighborhood of the boundary. 

Lemma 7.1. The kernel of Pi is smooth for x > and has a cornormal singularly at x = 0. 
Moreover, 

n 

p{x,y;y')^Y.'''udy,y') + 0{x^+^) 

1=0 

for some e > 0, where Ui G G"°°(aIT x dM; Hom{'f'^ , S'^{^T* M))) with singsupp(uj) G diag{dM) 
are distributions of finite order for all i. The asymptotic expansion of p holds in the sense that for 
any f G C°°{dM x dM; S^{°T*dM)), 

n 

Pi(/)=^u.(/) + 0(a;"+^), 

i=l 

with u^if) G C^{dM;S^{°T*M)). 

Since E{h,k) — E{Pi{hQ), P2{ko)) is a linear operator on Hq (g) fcg, we can denote its Schwartz 
kernel by 

e(m;y,y) G C^^{M x dM x dM;Hom{y^ ® r\ S'^(°r*M))). 

Lemma 7.2. The kernel of E is smooth in x > with only cornormal singularity at x — 0. 
Moreover, near the boundary 

n 

e{x, y; y,y)^Yl y, y) + 0(a;"+') 

i=l 

for some e > 0, where a G C-°° {{dMf; Hom.{y^ (g) , {°T*M))) with 

singsupp(e,) e {y ^ y} U {y ^ y} U {y ^ y} 

for all < i < n. The asymptotic expansion of e holds in the sense that for any /i,/2 G 
G°°((9M;ri), 

n 

E{fi ® /2) - ^ x'e^ifi /a) + 0(x"+^) 

i=0 

with Ciifi (g) /a) G G°°(5M; S'^{°T*dM)). 

Proof. First, for any e > and x > e, can choose a cut off function xi^) G C°°((Q,oo); [0, f]) such 
that x{^) = 1 for X > and = for x < ^. Then e{m;y,y) = E{x{x)p{m; y) , x{x)p{m; y)) is 
smooth since both xix)p{m; y) and y) are smooth, vanishing near the boundary, and E is 
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a differential oeprator with respect to each variable. The asymptotic expansion of e{x,y;y,y) can 
be compute directly. Since for any /i, /2 G C°°{dM\ 'V^), 



1=0 

with ui{fj) G C°°(9M; S'^{°T*dM)). The asymptotic property of E shows that 
E{Pih,Pif2) = E„+xE, + --- + x"E^ + 0{x^+') 

where E, e C°°{dM; S'^{°T*JI)) is a finite sum of quadratic forms on u;' (/2) for 1,1' < i 

and is a pseudodifferential operator with respect to each variable. So the kernel of Ei is a sum of 
distributions of the form 

/ _ w{y, y')ui {y', y)ui, {y' , y)dy' 
JdM 

with singsupp(ii;) e {y = y'}, singsupp(M() G {y' == y} and singsupp(u;/) e {y' = y}. li y ^_y, y ^ y 
and y y, there are small neighborhoods U oiy, U oi y and (7 of y, such that UOU — %,U\^U — % 
and C/ n J7 = 0. Then for y' in each of these small neighborhood, at most one of w,ui,uii can be 
singular. Hence the integral over such neighborhood will be well defined as a pairing of distribution 
with smooth function in the variable y' , smoothly depending on the parameter {y,y,y) G UUUU U . 
On dM — U — U — U, every term will be smooth and hence 

_ w{y,y')ui{y',y)ui'{y,y)dy' 

dM-U-U-U 

is smooth in {y,y,y) G UUUU U, too. □ 

Lemma 7.3. Suppose Ai(A) and A2(A) are two nonconstant analytic functions on C, satisfying 
Ai + A2 > n, and P G $0 ^'^''^'(Af) is a meromorphic family of pseudodifferential operators with 
kernel p^^p'^^F, where F & C°°(M XqM). Then for any k e R such that 3fiA2 - n + fc > and 
\i~ki No, and any u G C°°(M), 

P{x''u)^x''Gi^'X^^G2, 

where Gi,G2 G C°°(M), and 

GAdM = c(A)mU=o, G2\q-m = {x~^^P\x=o){u), 

where c(A) is some constant depending on A. 

IfX2-n + k^0, thenx'^Gi is replaced by x'' log xFi + x'^Gi with Fi,Gi G C°°(Af) and Fi =c'u 
for some constant d . 

// Ai - fc = / G No, then x^^G2 is replaced by x^^ \0gxF2 + x^^G2 with F2, G2 G C°°(M) and 

F2-^T,_,(a»(0,y) 

where Ti^j are differential operators of order I — j on the boundary for < j < I. 

Proof. Let ([/; y) be a small neighborhood in dM around p. Let (p{x) be a smooth cut off function 
such that for some e > 0, (l){x) = 1 if x < e and (j){x) = if a; > 2e. Let ip{y) be some smooth cut 
off funtion around p such that V'(y) — 1 if I?/ ~ pI < ^ ^-nd V'(y) = if I?/ ~ > 2e. Then 

P{x''u) = - <j))x''u) + P{<j){l ~ ^))x^u) + P{(i>^x^u). 

First, since (1 — (t))x^u G C^{M), clearly by Proposition 13.21 

P{{1 ~ (l))x''u){x,p) = 0{x^'). 

Secondly, 

Pm^ij)x''u){x,p)= r [ P,{p^^p'^-F{p,p',R,u;mx'){l-ij{y'))x"'u{x',y')dvolg. 

Jo J\y'-p\>e 
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Since j34p^^p'^^F{p,p',R,w))<j){x'){l - il){y')) e a;^ia;'^=C°°(M x M), clearly, 

P{4>{1 - ip)x''u){x,p) = 0{x^'). 

For the third part P{(j)tpx''u), we can choose / G C°°{dM), suppoted in {y G U : |y — p| < e}. Recall 
that z = y-y', R"^ = x'^ + x^ + uj = z/R, p = x/R and p' = x'/R. 

{P{,^i;x''u), f{y)dy) = {ttl. {p^^ p'^'^ FiTl{<p{x')w{y')x"'u{y'))), f{y)dy) 

= x'^ J J J R-^^-"+''p'^^--+'^F<f>{x')^{y')u{x',y')f{z + y')^dzdy'. 

Let v{x',y') = (j){x')tp{y')u{x' ,y'). Note that for any m G N, 

IlT=o{^d, - i)iRdn. + Xi+n-k- i)[R-^'-^+^ p'^^-^+'^FiR, p, p' ,u:,y')v(p' R,y')] 
= a;"'+ii?-^i-"+V'^^""+' (i^o + F,p'd,v{p'R, y') + ■ ■ ■ + F„p""5^^;(p'i?, y')), 
where Fj G C°°(M Xo M) for < i < m. Hence 

^"Loi^9^ - OM, + Ai - A; - i){x-^'P{cj)i:x''u), f{y)dy) 

= x^+^ j j f R->^--^+^p'>^^-^+\Y.Fip'^dlv{p'R,y'))f{z + y')%dzdy', 
for any m G N. This implies that for some m' independent of m, 



(P(0V^'=m), f{y)dy) = ^ x'Ga + x^' J2 "^^^ + 



^m+l+m' ^ 



i=0 



where Ga,Gi2 G C~(M) for < i < m. If fc < Ai, let 6 = z/x, t = x'/x and pp = [l + t^ + \e\^)-^ . 
Gio = ^^^J_l I P^/+^H^'-''+''F{xpp\pF,tpF,pFe,y')v{tx,y')f{xe + y')^dedy' 

= il f P^F'+'''t^'-''^''F{O,pp,tpF,pFe,y')^d0){v{O,y),f{y)dy). 

Note that pp is just the boundary defining function on the front face Fp with respect to its boundary 
surface B H Fp. So above integral is actually an integration over Fp which is well defined since 
Ai + A2 > n. If Ai > k, then letting i?o = -R|a;=o, 

G20 = J J J Ro^'+'p'^^-^+''F{R^,0,p',^,y')v{p'Ro,y')f{Rou> + y')^^dcody' 

= {{x-^'PU=o){x''v){y),f{y)dy) 

The two formulae can be extended meromorphically to all {A G C : A2(A) — n+fc > 0, Ai(A) + A2(A) > 
n} by analytic continuation. Adding the three parts together, then for m arbitrarily large, 

m m 

P{x^u) = x^^x'Gii + x^' ^x'Gi2 + 0(0;"'+^+'"') 

i=0 i=0 

where 

Gio = c{\)u{{),y), G20 = ix-^'Pl,=o){x''u){y): 

cW= I I p^p'^^H^'-^+''F{0,pF,tpF,pFe,y')^de. 

The distribution /Jq"^^^*^^^ on 1-dimentional space is a homogeneous distribution defined meromor- 
phically over Ai G C with poles of order 1 at Ai = A: + Nq. For / G No, the residue is 

Resx,=k+i{G2o) = J J p'^'-''+''d'j,^{F{Re,0,p',^,yMp'Ro,y-Roco))\R,=oydu; 

= 11 ^ p''-^-d'^^F{0,Q,p',u^,y')di^d'^u{Q,y))%d^ 

;2<i 
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where Ti^i^ are differential operators of order I — I2 on the boundary for < /2 < ^. Similarly 
^A2-n+fc-i pg^j-^ extended meromorphically to all A2 G C with poles of order 1 at A2 = n — fc — Z 
for Z e Nq. At A2 = n — fc, the residue is 

i?esA,=„_fc(GiO) = ^ p^^+^^i^(0, PF, 0, ppe, y')de. 

Herepi. = (l + l^p)-!. □ 

Let Hi : C°° {dM; S^{°T*M)) — > C°°{dM;'r^) be the projection map. 

Proposition 7.1. Given HQ,h(),ko as above. Then H = Pi(IIiHq) + Qg{hQ (g) fco) where Qg is a 
quadratic term with kernel q{m;y,y) G C"°°(M x dM x dM;Hom{'r^ y^,S^{"T*M))), defined 
by 

q{m;y,y)^ / Rg{m;m')e{m']y,y)dvolg{m'). 
Jm 

Moreover, in a neiborhood of the boundary, 

n 

<i{x, y, y,y) = Yl ^^^^(y; v, y) + 0(a;"+'=), 

1=0 

for some c> 0, where q^{y]y,y) e C-'=^{{dMf;Hom{'r^ ® , S^{°T* M)))) with 

singsupp(gi) c = y} U {y = y} U {y = y}. 

Proof For any 7 G ^""(Af; S^{°T*M)), let iJ' i?3(7). By SectionEl H' = a;"i/; + 0(x"+'=) for 
some 0, where H',, G C°°{dM; Y^). Then 

(i7,7)A/ = /" {HLgH' - LgHH')+E{h,k)H'dvolg 



HoHj^dvolgo + / E{h, k)Rg{j)dvolg 
Jm 

[ (Fi(nii/o) + k))jdvolg. 
Jm 



'M 

Hence Qg{ho (8) fco) = Rg{E{Piho, Pi fco)), whose kernel is 



'?("T.;y,y)= / Rg{m;m')e{m';y,y)dvolg{m'). 



gyiil, III, jcyiii , y, Llju,UKjLgy 

M 

Instead of dealing with the above integral directly, we introduce a parameter A G C as was done in 
Section [5] and deal with 

q{X,m;y,y) ^ / Rg{X,m;m')e{m';y,y)dvolg{m'). 
Jm 

for |A — n| < e with e > and small enough. 

Suppose Wi, Ui are small open sets in dM such that Wi C Wi C Ui for i = 1, 2, 3 and Ui, U2, U3 
are disjoint. Let U4 = dM — Uf^^W^i and {ipi}f^i is the partition of unity subject to the convering 
{Ui}j^i of dM. Let (j){x) be a smooth cut off function such that for some e > 0, 0(x) = 1 if a; < e 
and (j){x) = if x > 2e. Then for (y, y, y) G C/i x [/2 x f/s and < |A — rt| < e, 



q{X) = Rg{X){{l - 0)e) + ^i?3(A)(0V^,e). 



1=1 



First, by Lemma O (1 - 0)e G C;?°(Af x x 9M). Clearly, 

i?3(A)((l-0)e)=O(x^), 

smoothly depend on {y,y,y). Secondly, (j)ipie G C°°{M x x 9Af). By a proof similar to that of 
Lemma |7.3[ we have 

i?g(A)(0^ie) = Gil + x^Gi2 + 0(a;^+^), 
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where Gn, G12 E C°°{M x dM x dM) and c> 0. Moreover, 

Gi2U=o - {x-^Rg{\)\:c=o){Hie). 
Thirdly, for i = 2 or i = 3 or i = 4, Rg{X){x, y; x' , y')(t>{x')il;i{y') € .c/^^^-^«(M x M). Hence 

i?g(A)(0^A,e) = a;^G,2 + O(x^+") 
for some Gt2 €C°°(M x dM x 9M) and c> 0. Moreover, 

G.2U=o = e((x-^i?g(A)U=o)0V'»)- 
Let Gn = ELo 1^ + Oix"+^), and 

4 

qx = {x-^Rg{X)U=o){{l ~ (j))e + (l>^ie) + ^ e{{x~^Rg{X%^o)^ij,). 

1=2 

Then 

n 

(39) q(A) = ^g, + gA + 0(a;^+^) 

i=0 

for some c > 0. By the arbitrariness of Ui, 1/2,1/3, 

(40) singsupp((7i) U singsupp((7A) C {y ^ y} U {y ^ y} U {y ^ y}. 

According to Lemma [731 the results above can be extended meromorphically to {A G C : |A — n| < e} 
with a pole at A = n. The residue at A = n corresponds to the logx terms in q{n). However, this 
can not happen. By paring with ho{y) (g) ko{y) for ho.ko e C°°{dM; y^), 

Qg{ho^ko) = H-Pi{IliHo). 

Recall that the asymptotic expansion of neither H nor Pi{IliHo) contains any x" \ogx term, so the 
asymptotic expansion of Qg{hQ ® ko) can not contain any x" log a; term, either. Therefore (j39p and 
(|4n| hold for aU |A - n| < e. □ 

Corollary 7.1. gn'^^ = Tng^^'^'^ + ^n(3o*''^\ So^'^'') where Tn is a pseudodifferential operator of 
order n and Rn is quadratic form defined by kernel 

rn{y;y,y) € C-'^{{dMf\Horn{S^{T*dM)(g)S^{T*dM),S^{T*dM))) 

with singsupp(r„) C {y = y} U {?/ = U {y = y}. 

Finally, for the N point problem with iV > 3 and (3 — [1, ...,1) E N^, we can proceed in the same 
way as was for the 2-point problem, i.e. 

(1) Fix a suitable gauge such that f]g{d"g) = F{d'^g; 7 < a) for all a < /3 such that the linearized 
equation is linear and elliptic with repect to d^g; 

(2) Use boundary paring formula and the resovent of Lg to write out the kernel of nonliear 
Poisson operators and Scattering matrix. 

All the information about singularities of the Dirichlet-to-Neumann map we can get is from the 
analysis of the kernel in the second step. However, the first step is a key to make the problem 
solvable. 
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